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ABSTRACT 


Hugoniot  equation  of  state  data  in  the  pressure  range  5  to  250  kb 
have  been  obtained  for  quartzite,  sandstone,  calcite,  marble,  limestone, 
plagioclase,  and  basalt.  The  data  were  obtained  by  conventional 
shock  wave  techniques  which  measure  shock  and  associated  free-surface 
velocities.  Impedance  match  solutions  were  obtained  for  porous  rocks. 

High  values  of  the  Hugoniot  elastic  limit  were  observed  in  solid 
rocks — 50  to  100  kb  in  quartzite,  40  to  50  kb  in  feldspar  and  basalt, 
and  20  kb  in  calcite  and  marble.  Porous  rocks  show  considerably  reduced 
values  amounting  to  approximately  5  kb  in  sandstone  and  limestone. 

Four  phase  transitions,  indicated  by  multiple  shock  fronts,  are 
observed  in  calcite.  No  simple  relation  between  these  and  the  reported 
transitions  observed  in  static  experiments  is  evident.  An  indication  of 
a  phase  transition  in  sandstone  at  68  kb  was  also  obtained  and  is  inter¬ 
preted  as  the  quartz-stishovite  transition. 

Decay  of  the  elastic  precursor  wave  amplitude,  implying  a  time- 
dependent  yield  stress,  is  observed  in  quartzite,  calcite,  and  marble. 
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I  INTRODUCTION 


The  problem  of  distinguishing  an  underground  explosion  from  an 
earthquake  on  the  basis  of  seismic  recordings  is  a  complex  one  and  is 
being  studied  by  a  variety  of  approaches.  Since  the  two  phenomena  differ 
primarily  in  the  mode  of  energy  release  to  the  surrounding  rock,  it  is 
clear  that  a  complete  understanding  of  the  differences  in  seismic  behavior 
must  include  detailed  knowledge  of  shock  propagation  and  decay  in  the  vi¬ 
cinity  of  an  explosion. 

The  mathematical  description  of  shock  propagation  is  based,  in 
general,  on  three  relations,  t/iz.,  the  equation  of  continuity,  the  equa¬ 
tion  of  motion,  and  a  constitutive  relation  among  the  stress  and  strain 
tensors  and  the  internal  energy  of  the  medium.  The  first  two  are  uni¬ 
versal  and  independent  of  the  particular  medium  except  insofar  as  the 
properties  of  the  medium  may  dictate  whether  certain  simplifications  may 
reasonably  be  made.  The  third  must,  in  general,  be  determined  for  each 
medium  experimentally,  since  theory  is  not  yet  capable  of  providing 
adequate  predictions. 

In  this  report  we  describe  work  performed  during  the  past  year  to 
determine  the  constitutive  relation  for  several  simple  rocks  and  minerals 
in  the  pressure  range  0  to  250  kb.  These  are:  quartzite,  sandstone  and 
water- saturated  sandstone;  calcite,  marble,  and  limestone,  single  and 
polycrystalline  plagioclase,  and  basalt. 

Most  of  this  work  has  involved  measuring  the  llugoniot  equation  of 
state;  this  is  the  locus  of  equilibrium  atress-densi ty  states  attainable 
from  a  given  initial  state  through  a  single  shock  transition,  where 
“  stress”  refers  to  that  component  of  the  stress  tensor  normal  to  the 
ahock  front.  These  data  were  obtained  by  conventional  shock  wave  tech¬ 
niques  in  which  measurements  are  made  of  shock  wave  velocities  and 
associated  free-surface  velocities. 

An  important  feature  of  these  measurements  was  the  observation  of 
multiple  shock  fronts  due  to  yield  at  the  elastic  limit  for  each  material, 
and,  in  the  case  of  carbonate  rocks,  due  also  to  several  polymorphic  phase 
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transitions.  In  addition,  stress  relaxation  in  the  elastic  waves  was 
observed  in  quartzite,  marble  and  calcite;  some  indication  of  relaxa¬ 
tion  was  also  obtained  for  feldspar. 

A  major  goal  of  this  program  is  to  determine  the  constitutive  re¬ 
lation  for  simple  minerals  which  occur  commonly  in  the  earth's  crust. 

With  this  information  as  a  basis,  it  is  hoped  that  theoretical  models 
can  be  found  for  estimating  the  relations  obtaining  for  more  complicated 
rocks  composed  of  mixtures  of  theae  minerals.  The  work  thus  far  has 
dealt  primarily  with  minerals  and  mono-mineralic  rocks;  intensive  in¬ 
vestigation  of  mixtures  has  not  yet  been  undertaken.  However,  some 
preliminary  experiments  on  water- saturated  sandstone  were  performed. 

The  stress-density  relations  that  are  obtained  during  unloading  from  a 
shocked  state  are  more  difficult  to  measure  than  the  llugoniot  and,  while 
several  techniques  can,  in  principle,  provide  this  information,  these 
are  not  yet  well-established.  A  few  preliminary  experiments  for  this 
purpose  were  performed  on  sandstone  during  the  reporting  period.  These 
indicate  that  high  accuracy  is  required  in  order  to  obtain  meaningful 
results. 

Section  II  of  this  report  presents  a  summary  of  the  major  results. 

A  brief  description  of  experimental  techniques,  described  in  more  detail 
in  an  earlier  report,1  appears  in  Section  III.  Section  IV-A  deals  with 
the  experimental  results  and  interpretation  for  quartz  rocks.  Section  IV-B 
with  carbonate  rocks,  and  Section  IV-C  with  feldspars  and  basalt. 

The  results  of  static  compression  experiments  on  porous  rocks  are 
given  in  Section  V. 
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II  SUMMARY 


The  major  results  and  conclusions  may  be  listed  as: 

1.  Data  on  the  Rankine-Hugoniot  equation  of  state  have  been 
obtained  to  2S0  kb  for  quartzites  and  sandstones,  for 
calcite,  marble  and  limestone,  for  single  and  poly¬ 
crystalline  plagioclase,  and  for  basalt. 

2.  Quartz  rocks  have  an  elastic  wave  amplitude  which  ranges  from 
5  kb  in  sandstone  to  40  kb  in  Sioux  quartzite,  to  60  kb  in 
Eureka  quartzite,  and  to  100  kb  novaculite.  The  elastic  wave 
amplitude  decreases  as  the  porosity  increases.  In  addition, 
decay  of  the  elastic  wave  amplitude  with  distance  of  travel  is 
observed  in  the  quartzites.  The  temperatures  attained  in 
porous  sandstones,*  when  shock  loaded,  are  much  higher  than 
those  attained  in  quartz  or  quartzites.  Temperatures  of 
900°K  are  estimated  in  Coconino  sandstone,  shocked  to 

95  kb.  A  possible  interpretation  for  some  of  the  results 
for  sandstone  is  that  the  quartz-stiahovite  transition  is 
observed. 

3.  Calcite  and  carbonate  rocks  have  an  elastic  wave  amplitude 
that  varies  from  about  20  kb  in  calcite,  to  about  12  kb  in 
Yule  and  Vermont  marble,  to  about  10  kb  in  Solenhofen  lime¬ 
stone,  and  to  about  5  kb  in  Spergen  limestone.  Besides  the 
Hugoniot  elastic  limit,  there  are  three  additional  transi¬ 
tions  in  calcite  which  are  indicated  by  multiple  shock 
fronts.  These  occur  at  about  30,  45,  and  95  kb;  values  at 
each  transition  vary  slightly  depending  upon  crystal 
orientation.  The  45  and  95  kb  transitions  also  occur  in 
marble;  no  evidence  was  obtained  for  the  30  kb  transition. 
However,  a  22  kb  transition  is  observed  in  marble  (the 
Hugoniot  elastic  limit  is  12  kb).  Calcite  and  Vermont 
marble  show  decay  of  the  elastic  wave  amplitude  with  dis¬ 
tance  of  travel. 

4.  Plagioclase  and  basalt  have  Hugoniot  elastic  limits  of 
about  40  kb  and  have  no  other  observed  transitions  for 
pressures  as  high  as  250  kb.  There  is  some  indication 
that  the  elastic  amplitude  decays  with  distance  of  travel. 

5.  Transitions  recorded  by  Adadurov,  Balashov,  and  Dremin2  in 
marble  and  by  Hughes  and  McQueen1  in  gabbro  may  require 
some  reinterpretation.  Similar  transitions  in  calcite  and 
basalt  were  not  observed. 


*  •■<»>«•  hpiiMi  obtoiaed  will  tbe  nl|i  (h{,  1,  Pig*.  IS  mi  1*,  uS 

Ttkl*  IV,  pagao  ^  sad  W  «•  »°w  kaowa  to  bo  incorrect  above  IS  bb. 
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a.  The  transition  at  ISO  kb  in  marble  reported  by 
Adadurov  et  al.  is  probably  the  transition  we 
recorded  at  about  95  kb,  150  kb  being  the  upper 
limit  at  which  the  multiple  shock  from  the  95  kb 
transition  may  be  observed. 

b.  The  estimated  transition  at  150  kb  in  gabbro 
reported  by  Hughes  and  McQueen  may  be  the  re¬ 
sult  of  elastic  yield  at  about  50  kb.  This  is 
indicated  by  data  using  basalt— a  rock  similar 
to  gabbro,  and  plagioclase  —  a  mineral  occurring 
in  sizable  amounts  in  gabbro. 

6.  Yield  stress  in  static  one-dimensional  strain  experiments 
for  sandstone  is  nearly  the  same  value  as  that  obtained 
in  shock  experiments.  This  correlation  is  not  observed 
for  limestones. 
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Ill  SHOCK -WAVE  EXPERIMENTS 


The  methods  used  in  determining  Hugoniot  equations  of  state  are 
discussed  in  Semiannual  Technical  Report  No.  2,  April  15,  1962.  These 
methods,  as  well  as  methods  reviewed  by  Rice  et  al . ,  are  baaed  on  the 
Rankine-Hugoniot  jump  conditions  expressing  conservation  of  mass,  momentum, 
and  energy  across  a  shock  front. 


(1) 
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where  V  is  specific  volume,  u  is  mass  velocity,  U  is  shock  velocity,  cr  is 
compressive  stress  normal  to  the  shock  front,  and  E  is  specific  internsl 
energy.  All  shock  front  and  particle  velocities  used  in  this  report 
refer  to  laboratory  coordinates.  Subscripts  "0”  refer  to  the  state  ahead 
of  the  shock  front,  and  subscripts  “  1  ”  to  the  state  immediately  behind 
the  front. 

Assuming  the  state  ahead  of  a  shock  front  to  be  known,  there  are 
five  unknown  quantities  and  three  relations  to  describe  the  state  behind 
the  shock.  Hence,  measurement  of  any  v.wo  allows  all  quantities  to  be 
computed.  In  the  method  used,  the  shock  wave  velocity  and  the  velocity 
from  a  free  surface  are  measured.  These,  together  with  a  relation  between 
free-surface  velocity  and  particle  velocity  prior  to  reflection,  are  suf¬ 
ficient  to  determine  the  stress- density-energy  state  behind  the  shock 
front.  A  series  of  measurements  with  different  shock  pressures  determines 
the  Hugoniot  equation  of  state. 

The  particle  velocity  prior  to  reflection  is,  to  an  adequate  approxi¬ 
mation,  one-half  the  free  surface  velocity.  This  approximation  is  bssed 
on  the  condition  that  entropy  changes  are  small.1  In  essence,  this  means 
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that  the  unloading  atatea  ahould  follow  the  mirror  image  of  the  Hugoniot 
atatea  represented  in  the  pressure-particle  velocity  plane.  This  assump¬ 
tion  is  satisfactory  for  crystals  and  non-porous  rock,  but  is  not  for 
porous  rocks,  i.e-j  limestones  and  sandstones. 

The  pressure  region  from  a  few  kilobars  to  a  few  hundred  kilobars  is 
the  region  in  which  the  elastic  strength  and  phase  transitions  affect  the 
shock.  This  is  the  region  in  which  we  have  made  our  measurements.  Such 
measurements  in  rocks  are  difficult  because  of  the  wide  range  in  ampli¬ 
tudes  of  elastic  precursors,  and  because  rock  surfaces  are  nonconducting, 
brittle,  and  poor  reflectors.  Occurrence  of  phase  transitions  makes 
problems  more  difficult.  Desirable  techniques  are  those  that  measure 
continuously  ahock  and  free-surface  velocity.  The  techniques  that  were 
used  are  briefly  described  below. 

A.  INCLINED  MIRROR  METHOD 

The  inclined  mirror  method  was  used  in  moat  experiments.  More  detail 
is  presented  in  Semiannual  Technical  Report  No.  2,  April  15,  1962,  and  in 
a  paper  by  G.  R.  Fowles. * 

The  rock  specimens  are  usually  about  54-inch  or  54-inch  thick  by 
154-inches  square.  A  specimen  is  cemented  to  a  backing  plate  and  a  plane 
shock  wave  is  driven  through  the  backing  plate  and  into  the  specimen.  An 
explosive  assembly,*  which  produces  a  plane  detonation  wave,  is  cemented 
to  one  surface  of  the  backing  plate;  specimen  and  mirrors  are  cemented  to 
the  other  surface  (Fig.  1  and  2).  The  experiment  is  designed  to  determine 
shock  and  particle  velocity.  Shock  wave  velocity  is  measured  by  observing 
the  transit  time  of  a  shock  wave  through  a  known  thickness  of  sample. 
Particle  velocity  ia  determined  by  measuring  the  free-surface  velocity 
produced  upon  interaction  of  a  ahock  wave  with  a  free  boundary  and  con¬ 
verting  to  particle  velocity  as  discussed  above. 

Figure  2  is  a  simplified  diagram  of  the  interaction  of  the  assembly 
with  a  plane  shock  wave.  The  shock  front,  Sj  in  the  backing  plate,  gener¬ 
ated  by  the  explosive  lens  travels  with  a  velocity  U j.  Mirrors  along  with 
the  aluminized  lower  surface  of  the  specimen  (if  the  specimen  is  trans¬ 
parent)  record  the  arrival  of  the  shock  wave  at  the  upper  surface  of  the 


* 

Maaafactarad  by  Mato*  and  Haagar,  lac.  (  Aaiarillo,  Taaaa. 
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FIG.  1  INCLINED  MIRROR  ASSEMBLY 

(Vacuum  chamber  not  shown) 


becking  plete.  This  is  recorded  as  a  change  in  the  intensity  of  reflected 
light  when  the  shock  wave  destroys  the  alumini  s**  l  surface  of  the  mirrors. 

When  the  shock  wave  traveling  in  the  backing  plate  reaches  the 
boundary  between  the  plate  and  specimen,  it  interacts  to  transmit  a  shock 
wave  of  different  pressure  into  the  specimen.  Assume  the  shock  wave  in 
the  specimen  is  unstable  and  forms  a  double  shock,  and  $j,  indicating 
a  phase  change  or  failure  at  the  elastic  limit.  The  two  shock  waves, 
and  in  the  sample  travel  with  velocities  t! f  and  f/j.  The  reflection 
of  the  shock  wave  with  velocity  f/j  at  the  surface  of  the  sample  acceler¬ 
ates  the  surface  to  velocity  u^j,  which  is  much  slower  than  the  shock 
wave  velocity;  then  the  shock  wave,  with  velocity  f.'j,  overtakes  the 
free  surface  and  further  accelerates  it  to  velocity  u^j. 
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FIG.  2 


Initially,  the  inclined  mirror  on  the  sample  is  turned  by  the  impact 
of  the  surface  traveling  with  velocity  u^j,  and  if  the  angle  between 
mirror  and  specimen  surface  is  a,  the  point  of  collision  travels  along 
the  mirror  with  apparent  velocity  U •  u^/ain  <*.  When  the  second  shock 
accelerates  the  free  surface  to  velocity  u^j,  the  apparent  velocity,  Uai, 
is  increased  to  u^j/sin  tt.  These  apparent  velocities  are  recorded  with 
a  streak  camera.  Thus  distances  indicated  by  t  on  the  film  in  Fig.  2 
are  related  to  shock  velocities.  The  angles,  y,  are  related  to  particle 
velocities.  From  the  measured  shock  and  particle  velocities  the  other 
quantities  describing  the  shocked  state  can  be  computed.  Note,  too, 
that  this  technique  provides  a  continuous  measurement  of  free-surface 
velocity. 

B.  IMPEDANCE  MATCH  SOLUTION 

The  impedance-match  solution  for  the  stress-particle  velocity  state 
in  the  specimen  is  obtained  by  the  method  shown  in  Fig.  3.  The  curve  OA 

is  the  previously  determined  locus  of  points  which  can  be  attained  by 

passing  a  shock  through  the  backing  plate.  The  point  A  represents  the 
state  behind  the  shock  which  is  incident  on  the  boundary  between  the 
backing  plate  and  the  specimen.  The  particle  velocity,  u,  is  obtained 
from  the  free-surface  velocity  of  the  backing  plate  and  is  measured  with 
an  inclined  mirror.  A  reflected  wave  is  generated  when  the  shock  reaches 
the  interface  between  the  backing  plate  and  specimen.  The  states  which 
can  be  attained  by  the  reflected  wave  can  be  calculated  from  the  known 
Hugoniot  of  the  backing  plate;  they  are  very  close  to  AB,  a  mirror  image 

of  the  curve  OA.  Since  normal  stress  snd  particle  velocity  must  be  con¬ 

tinuous  across  the  interface,  the  state  in  the  sample  must  lie  on  AB  . 

From  the  equation  for  conservation  of  momentum,  *pQU1ul.  the  state  in 
the  specimen  must  also  be  on  a  straight  line  with  a  slope  ( c^/u,  )  •  p#(/, . 
The  intersection  of  this  line  with  AB  determines  B,  the  state  in  the 
specimen.  The  associated  compression  and  specific  energy  may  be  calcu¬ 
lated  from  the  equations  for  conservation  of  mass  and  energy  across  the 
shock  front. 

Sometimes  free-surface  velocities  from  two  different  metals  placed 
on  top  of  the  specimen,  or  the  free-surface  velocity  from  the  backing 
plate  and  the  particle  velocity  from  a  metal  on  top  of  the  specimen  are 
measured.  The  particle  velocities  are  measured  with  inclined  mirrors, 
optical  lever  techniques}  or  quartz  transducers.  For  example,  assume  the 
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FIG.  3  IMPEDANCE  MATCH  SOLUTION 


shock  velocity  in  the  sample  and  particle  velocities  of  two  metals  on  top 
of  the  specimen  have  been  measured.  The  f ree-sur face  approximation  is 
used  to  obtain  particle  velocities  in  the  two  metal  mirrors  on  the  rock 
specimen,  giving  points  A  and  B  on  the  known  a  -  u  loci  of  these  metals 
(Fig.  4).  A  lies  on  the  locus  of  states  attainable  by  reflected  shocks 
in  the  rock  and  B  lies  on  the  locus  of  reflected  rarefactions  from  the 
rock.  If  A  and  B  are  close  to  one  another,  the  straight  line  AB  is  a 
good  approximation  to  the  cross  curve  of  reflected  wave  states  of  the 
rock  Hugoniot.  The  intersection  at  C  of  AB  with  the  line  from  the  origin 
of  alope  p0U  determines  the  state  from  the  rock. 

C.  QUARTZ  TRANSDUCER  METHOD 

The  measurement  of  shock  wave  profiles  in  aluminum,  various  steels, 
iron,  and  barium  titanate  by  means  of  £-cut  quartz  transducers  has  been 
described  by  Jones  et  al.,1  and  Halpin  et  al ,*  These  quartz  transducers 
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FIG.  4  IMPEDANCE  MATCH  SOLUTION 


operate  up  to  a  stresa  level  of  21  kb.  A  longitudinal  atreaa  difference 
parallel  to  the  X  direction  generatea  a  displacement  current  by  piezo¬ 
electric  polarisation.  A  longitudinal  atreaa  wave  impinging  on  the 
quarts  transducer  (Fig.  5)  generates  a  current.' 

Akv 

»(*)  •  ItrAO.t)  -cr(h.l)]  (♦) 

A. 

where  A  is  the  effective  area  of  the  transducer,  ii  the  transducer 
thickness,  vf  the  longitudinal  elastic-wave  velocity  in  quarts  in  the  X 
direction,  the  longitudinal  stress  in  the  quarts  (a  function  of  position 
and  time),  and  k  ia  the  piesoelectric  coefficient  for  strain  in  the  X  di¬ 
rection.  By  using  Eq.  (4)  it  is  assumed  that  the  strain  and  the  displace¬ 
ment  field  are  one-dinensional ,  and  the  dielectric  constant,  piesoelectric 
coefficient,  and  longitudinal  wave  velocity  are  independent  of  the  stress. 
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FIG.  5  EXPERIMENTAL  CONFIGURATION  FOR  QUARTZ  TRANSDUCER 
STRESS-TIME  PROFILES  IN  POROUS  ROCK 


The  stress-wave  generator,  specimen,  and  quartz  transducer  dimensions 
are  adjusted  so  that  the  duration  of  the  stress-wave  profile  of  interest 
is  less  than  the  one-way  travel  time  of  a  stress  wave  through  the  quartz 
tranaducer.  For  this  arrangement  the  time  interval  of  interest  is, 

0  <  t  <  h/v  ,  and 

9  9 

the  atreas  on  the  side  of  the  quartz  away  from  the  apecimen, 

crf(/i,  t)  ■  O 

Then  Eq.  (4)  may  be  written  in  the  form 

h  i(t) 

or. AO,  t)  -  -f- -  (5) 

*  On 


where  the  current  generated  by  the  quartz  transducer  is  now  proportional 
to  the  stress  at  the  specimen-quartz  interface. 

In  practice  the  one-dimensional  geometry  implicitly  required  in  the 
above  discussion  is  experimentally  obtained  by  using  only  the  central 
portion  of  the  quarts  tranaducer.  A  guard- ring  arrangement  in  which  the 
electrode  covers  only  the  central  region  of  the  disk  was  used  in  one  of 
the  present  experiments  (Section  IV-A-2).  Specific  electrode  arrange¬ 
ments  are  discussed  by  llalpin  «t  ol.* 
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If,  as  depicted  in  Fig.  S,  the  entire  surface  of  the  quartz  ia  used 
to  obtain  the  current  from  the  transducer,  significant  errors  due  to 
edge  effects  are  introduced.  In  the  present  experiment  these  errors  may 
be  evaluated  from  the  work  of  Jones  et  alJ  For  a  transducer  diameter-to- 
thickness  iatio  of  5  to  1,  and  in  the  recording  time  interval  of  0.00  to 
0.35  fiaec,  the  stress  indicated  by  the  quartz  transducer  using  Eq.  (5)  is 
approximately  3%  too  low. 

The  stress-time  record  in  the  quartz  transducer  requires  some  ad¬ 
ditional  assumptions  if  we  wish  to  obtain  stress  as  a  function  of  time 
for  the  specimen.  Consider  a  point  A  on  the  pressure-particle  velocity 
curve  of  JY-cut  quartz  [Fig.  6(a)]  where  A  represents  the  state  in  the 
quartz  at  the  quartz- specimen  interface  at  a  given  instant  of  time. 

In  order  that  the  state  A  in  the  quartz  be  related  to  the  state  B 
in  the  specimen  it  ia  necessary  that  the  shock  reflection  path  AB  be 
known.  Generally  the  shock  leflection  curve  as  well  as  the  shock 
compression  curve  OB  are  unknown.  For  a  specimen  of  higher  shock  imped¬ 
ance  (not  shown  in  Fig.  6)  than  the  quartz,  a  rarefaction  wave  will  con¬ 
nect  the  state  in  the  specimen  to  the  state  in  the  quartz.  For  a  multiple  step 
shock  front,  the  impedance  match  solution  is  computed  stepwise.  Through  the 
use  of  the  appropriate  electronic  timing,  the  relative  shock  velocity  (I/~u£), 


(o)  <b) 

FIG.  6  GRAPHICAL  SOLUTION,  QUARTZ  TRANSDUCER  DATA 

(a)  Erroneous  result  obtained  at  Br;  BA  is  actuol  shock 
reflection  path  rather  than  mirror  approximation  B' A; 

(b)  State  A'  is  observed  by  transducer  mounted  on  material 
of  stress-particle  velocity  curve  OM;  C  on  OM  corresponds 
to  A*;  for  solution,  ACB  is  drawn  as  approximation  to 
shock  reflection  curve  AB 
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appropriate  to  the  stress  level  A  and,  hence  Section  III-B,  the  slope  of 
the  line  through  <rg  is  known. 

In  the  present  experiments,  on  dry  and  wet  sandstone,  shock  reflection 
such  as  depicted  in  Fig.  6(a),  occurs  since  the  sandstones  are  of  lower 
shock  impedance  than  the  X-cut  quartz.  For  an  initially  porous  material 
such  as  sandstone,  the  use  of  the  mirror  image  line,  AB'  [Fig.  6(a)],  can 
give  a  poor  approximation  of  the  specimen  shock  reflection  curve.  One 
technique  for  obtaining  the  shock  reflection  curve  AB  and  the  state  B  in 
the  specimen  is  to  place  another  material  of  known  equation  of  state  on  top 
of  the  specimen  (Section  II I • B) .  Another  pressure  versus  time  record  is 
obtained  using  a  second  quartz  transducer  placed  on  the  material  surface  as 
depicted  in  Fig.  5.  At  a  given  time,  the  state  A'  is  obtained  in  the  quartz 
transducer  from  the  material  of  known  equation  of  state.  In  the  present 
experiment,  magnesium  is  used  for  the  material  of  known  equation  of  state. 
The  mirror  image  of  the  magnesium  pressure-particle  velocity  loading  curve 
is  uaed  as  an  approximation  of  the  shock  reflection  curve.  The  state  C  in 
the  magnesium  corresponding  to  the  state  A'  in  the  quartz  may  be  obtained 
graphically  [Fig.  6(b)],  In  order  to  relate  the  corresponding  states  in  the 
stress-time  profile  observed  from  the  two  transducers,  it  is  convenient  to 
assume  that  the  atress  profile  does  not  change  shape  in  the  magnesium,  ».e., 

t0  +  ~J  (6) 

Here  cr m  is  the  stress  in  the  megnesium,  a  function  of  position  xa  and 
time  t,  and  vm  is  the  magnesium  shock  velocity  which  is  assumed  constant.* 
The  points  A  and  A',  which  correspond  to  the  same  state  B  in  the  specimen, 
can  be  correlated  on  the  basis  of  equal  time,  as  measured  from  the  begin¬ 
ning  of  the  shock-wave  profile.  By  assuming  that  the  appropriate  segment 
of  the  reflected  shock  pressure-particle  velocity  curve  in  the  specimen 
may  be  approximated  by  a  straight  line,  the  line  AC  is  constructed.  The 
intersection  of  AC  with  the  line  of  slope  pt(U  ~  uf),  where  pf  and  uf 
are  the  density  and  particle  velocity  at  the  llugoniot  elastic  limit,  de¬ 
termines  [Fig.  6(b)]  the  state  in  the  specimen. 


*•*•»  10  kfc  tk*  ckaaga  ia  akack  aalaaity  iaaa|aa»iaa  witk  atrasa  ia  aa  tka  arSar  at  O.M/kk. 
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IV  RESULTS 


A.  QUARTZ  ROCKS 

1.  Quartz  and  Quartzites 

The  agreement  between  Fowles* 6  and  Wackerle’ a*  ahock  compression 
data  for  crystal  quartz  demonstrates  that  the  Hugoniot  is  repeatable  within 
experimental  accuracy.  The  Hugoniot  of  novaculite  and  quartzite  is  also 
reproducible  and  in  general  agreement  with  average  values  for  crystal 
quartz  (Fig.  7). 

For  different  cuts  of  quartz  it  is  not  clear  how  Hugoniot  data  are 
beat  represented  in  the  stress- volume  or  stress-particle  velocity  plane  so 
that  the  results  may  be  compared  to  those  for  polycrystalline  quarts.  Com¬ 
parison  with  Bridgman* a1#  di latometric  isothermal  measurements  is  also  use¬ 
ful  in  interpretation  of  shock  wave  results. 

Wackerle  has  calculated  the  shock  and  particle  velocity  in  quartz  ap¬ 
propriate  for  hypothetical  shock  transition  directly  from  zero  stress  to 
states  normally  achieved  only  via  a  second  wave.  These  hypothetical  shock 
and  particle  velocities  are  combined  for  X- ,  Y- ,  and  7-cuts  and  the  result 
is  analytically  fit  to  straight  lines  in  different  stress  ranges.  Below 
144  kb,  Wackerle* a  straight  line  was  fitted  to  a  Murnaghan  equation  of 
state  of  the  form  of  Eq.  (16),  the  result  is  shown  in  Fig.  7.  The  Bridgman 
data  has  also  been  put  in  the  form  of  Eq.  (16)  and  is  shown  in  Fig.  7.  The 
agreement  among  the  reduced  Hugoniot  data  of  Wackerle,  the  static  Bridgman 
data,  and  the  states  achieved  by  the  second  wave  in  polycrystalline  quartz, 
here  represented  by  Arkansas  novaculite,  is  noteworthy. 

Although  the  general  behavior  of  single  and  polycrystsl line  nonporous 
quarts  rock  at  relatively  low  pressures  is  thought  to  be  understood, 
several  phenomena  at  high  stresses  have  not  heen  explained.  These  include: 

(1)  The  differences  in  Hugoniots  above  the  elastic  limit  for 
various  crystal  directions  of  quartz  cannot  be  explained 
in  terms  of  present  models  of  material  behavior  under 
shock  compression. 
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(2)  The  question  of  shock  induced  polymorphism  in  quartz, 
which  is  suggested  by  Wsckerle  and  recently  discussed 
by  McQueen  <t  oi.u 

(3)  The  stress  relaxation  phenomenon  and  the  variation  of 
the  precursor  wave  amplitude  in  different  polycrystal¬ 
line  quartz  materials  require  further  investigation. 

2.  Sandstones 

The  shock  compression  of  Coconino  sandstone  has  been  examined  up  to 
186  kb.  The  results  for  the  Coconino  sandstone  and  those  from  several 


FIG.  7  STRESS-VOLUME  RELATIONS,  NONPOROUS  QUARTZ 
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experiment*  below  25  kb  using  Massillon  and  St.  Peter  sandstone  sre 
shown  in  Figs.  8  and  9. 

Coconino  sandstone  is  composed  of  97%  quartz  grains  ranging  in 
diameter  from  0.12  to  0.15  mm.  The  impurities,  3%  of  the  solid  volume, 
are  principally  microcline  and  orthoclase.  Porosity  of  this  rock  varies 
from  22.5-  to  24.1%.  The  present  data  (Table  I)  have  been  obtained  from 
material  in  which  the  beat  estimate  of  the  porosity  is  24.1%. 

Massillon  sandstone  is  composed  of  94%  (of  the  solid  volume)  quartz 
grains  ranging  from  0.01  to  0.5  mm.  The  6%  impurities  are  chiefly  feld¬ 
spar,  ircn  oxide,  and  clay.  Variations  in  color  suggest  the  clay  and 
iron  oxide  content  changes  in  the  rock.  The  rock  has  a  porosity  of  23.5%. 

St.  Peter  sandstone  is  composed  of  99%  (of  the  solid  volume)  of 
rounded  grains  ranging  in  diameter  from  0.05  to  0.4  mm.  Impurities  con¬ 
sisting  of  feldspar,  sphene,  zircon,  tourmaline,  and  clinozoisite  amount 
to  less  than  0.5%.  The  porosity  is  19.5%. 

The  Hugoniots  of  Coconino  and  Massillon  sandstones  do  not  appear  to 
differ  significantly  for  shocks  below  25  kb.  The  greater  shock  impedance 
of  the  St.  Peter  sandstone  is  probably  due  to  the  lower  porosity.  At 
zero  pressure  the  longitudinal  wave  velocity  determined  by  conventional 
ultrasonic  travel-time  methods  is  3.254  ±  0.178  mm//xsec,  3.183  i  0.216  mm/ 
Msec,  and  3.333  i  0.103  mm/^aec,  for  Massillon,  St.  Peter,  and  Coconino 
sandstone,  respectively.  These  velocities  are  in  fair  agreement  with 
those  of  the  elastic  precursor  listed  in  Table  I.  The  relatively  lower 
porosity  of  St.  Peter  sandstone  is  thought  to  produce  the  relatively  high 
precursor  velocity  of  approximately  3.63  mm/M*ec;  however,  this  high  value 
has  not  been  observed  in  the  longitudinal  ultrasonic  velocities. 

For  shock  states  up  to  the  Hugoniot  elastic  limit  (Table  I,  Figs.  8 
and  9)  the  inclined  mirror  technique  was  used  to  obtain  free-surface 
velocity.  The  free-surface  approximation  was  applied  in  order  that  the 
initial  state  for  the  following  second  shock  could  be  calculated.  In  the 
experiments  in  which  two  waves  were  observed,  the  final  state  was  calcu¬ 
lated  using  the  impedance-match  solution  with  the  initial  state  assumed  to 
be  that  of  the  precursor.  The  free-surface  velocities  which  correspond  to 
the  final  state,  for  shocks  in  which  two  waves  were  present,  are  indicated 
in  Fig.  8.  These  data  were  obtained  from  the  slope  of  the  inclined  mirror 
cut-off  on  the  streak  camera  records;  the  records  of  Figs.  10  and  11  being 
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FIG.  8  HUGONIOT,  SANDSTONE 


FIG.  9  STRESS-PARTICLE  VELOCITY,  SANDSTONE 


FIRST  FRONT 


SPECIMEN  TTPE 

SPECIMEN 

THICKNESS 

(aa) 

54  (Fraa- 

Sur  fact 
Valocity 
(aa/paac) 

Shock 

Valocity 

(aa/paac) 

Streii 

(kb) 

Voluae 

(cc/*) 

Coconino  Sandstone 

6.32 

6.35 

Coconino  Sandatone 

6.63 

6.43 

Coconino  Sandatone 

6.35 

6.32 

Coconino  Sandatone 

5.07 

0.156 

2.622 

(8.0) 

0.480 

9.29 

0.074 

2.853 

4.1 

0.497 

Maaaillon  Sandatone 

6.60 

0.109 

2.861 

(6.4) 

0.472 

Coconino  Sandstone 

6.61 

12.79 

o.oot 

Coconino  Sandatone 

6.34 

0.151 

2.705 

(8) 

0.453 

9.58 

0.086 

3.027 

5.1 

0.495 

Maaaillon  Sandstone 

5.21 

9.59 

0.080 

0.120 

2.540 

2.357 

4.2 

(6.9) 

0.475 

0.4697 

St.  Peter  Sandatone 

6.35 

9.57 

o.oot 

0.142 

0.058 

3.579 

3.670 

(10.9) 

4.5 

0.449 

0.460 

Coconino  Sandatone 

6.35 
9.53. 
0.00  1 

St.  Peter  Sandstone 
(Water  Saturated) 

7.23 

4.81. 

o.oot 

Coconino  Sandatone 
(Water  Saturated) 

7.20 

4.81 

o.oot 

Coconino  Sandatone 
(Water  Saturated) 

9.54 

0.021 

Paiticle 

Velocity 

3.222 

1.5  t  0.2 

0.447 

Coconino  Sandstone 

8.37 

0.043 

Particle 

Velocity 

1.881 

1.8  i  0.2 

0.498 

Coconino  Sandatone 

6.38 

9.58 

St.  Peter  Sandatone 

6.49 

12.69 

NOTE;  p#  tram  9197,  9106,  9209  aaad  for  all  cilnluiMi. 
(  )  iilieitM  tm  riliitli  data. 
tLiatar  Extrapolation. 

Snparatt  obaaraatioa. 


Table  I 

HUG0N10T  DATA  FOR  SANDSTONE 


FIRST  FRONT 

FINAL  STATE 

DRIVER 

N  TYPE 

SPECIMEN 

THICKNESS 

H  (Free- 

Shock 

Str«u 

Vol UNO 

Particle 
Velocity 
(aa/^aec  ) 

Free- 

Surface 

Shock 

Streaa 

Vol use 

EXPLOSIVE  SYSTEM 
P-60  ♦ 

Material 

Free- 

Surface 

(an) 

Velocity 
(»»//!»•  c) 

Velocity 

(aa/M*«c) 

<kb) 

(«*/») 

1  Velocity 
(*»//i.sec ) 

c  ) 

(kb) 

(ce/g) 

( inch ) 

Velocity 
(a»/|isec ) 

andatone 

6.32 

-in 

2.50 

4.126 

144 

0.289 

1"  Comp  B 

0.5  A1 

2.568 

6.35 

xm 

2.24 

4.039 

124 

0.311 

1"  Comp  B 

1  A1 

2.206 

in da tone 

6.63 

1.26 

1.95 

3.285 

83 

0.309 

2"  Baratol 

0.25  A1 

1.77 

6.43 

1.11 

1.63 

3.141 

70 

0.323 

2"  Baratol 

0.75  A1 

1.54 

■ndatone 

6.35 

1.09 

1.88 

3.321 

73 

0.337 

0.  5  Lucite 

2.46 

6.32 

1.02 

1.99 

3.136 

64 

0.337 

1  Lucite 

2.25 

indatone 

5.07 

0.156 

2.622 

(8.0) 

0.480 

0.371 

0.699 

2.305 

18 

0.402 

P-80  (instead 

0.5 

Brass  + 

0.836 

.  of  P-60) 

1 

Lucite 

9.29 

0.074 

2.853 

4.1 

0.497 

0.371 

0.499 

2.385 

18 

0.402 

1 

1 

Sandstone 

6.60 

0.109 

2.861 

(6.4) 

0.472 

0.354 

(0.732) 

2.553 

19 

0.406 

♦ 

i 

< 

1 

indatone 

6.61 

1 

4.481 
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Table  I 

IIUGONIOT  DATA  FOR  SANDSTONE 


FINAL  STATE 

DRIVED 

article 

al oci ty 
m/fxttc ) 

Free- 

Surface 

Shock 
Velocity 
(am/^aec  ) 

Street 

Voluee 

EXPLOSIVE  SYSTEM 
P-60  ♦ 

Material 

Free* 

Surface 

SHOT  NO. 

NOTES 

Velocity 

(ee/uaac) 

(4b) 

(cc/f ) 

( inch ) 

Ve  loci ty 
imm/fine ) 

1.74 

4.126 

' 

144 

1"  Comp  B 

0.5  A1 

2.568 

8974 

1.53 

4.039 

124 

0.311 

1"  Comp  B 

1  A1 
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8974 

1.26 

3.285 

83 

0.309 

2*  Raratol 

0.25  A1 

1.77 

8975 
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70 

0.323 
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0.75  A1 

1.54 
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0.337 
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64 
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18 
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9107 

i  of  P-60) 

1 
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0.499 

2.385 

18 

0.402 

1 
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2.553 

19 

0.406 

♦ 
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FIG.  10  INCLINED  MIRROR  CUTOFF 
COCONINO  SANDSTONE, 
DOUBLE  WAVE 
First  wave  4.1  kb;  second 
wave  18  kb 


FIG.  11  INCLINED  MIRROR  CUTOFF  COCONINO  SANDSTONE, 
SINGLE  WAVE 

Single  wave,  showing  lateral  relaxation  wave  interaction 
width,  W 
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typical.  The  general  indiatinctneas  and  unevenneaa  in  Fig.  10  are  thought 
to  produce  some  of  the  scatter  in  the  data,  particularly  at  free-surface 
velocities  of  less  than  2  mm/juaec.  Shock  attenuation  and  edge  effects  in 
porous  rocks  have  limited  most  of  the  experiments  to  specimens  between 
S-  and  13-mm  thick  with  a  6-  or  8-inch  diameter  explosive  lens.  For  these 
dimensions  the  grain  size  of  the  rocks,  particularly  when  on  the  order  of 
0.5  mm,  undoubtedly  affects  the  reproducibility  of  the  free-surface  veloci  ty 
measurements. 

The  porous  sandstones  studied  may  be  characterized  as  displaying  a 
Hugoniot  elastic  limit  in  the  vicinity  of  5  kb.*  Above  the  elastic  limit 
and  up  to  25  kb,  successfully  higher  stress  states  are  achieved  by  irre¬ 
versible  crushing  of  the  sandstone  until  the  quartz  density  is  reached.* 
Although  the  free  surface  velocity  data  are  not  entirely  self-consistent 
(Fig.  8,  see  also  Section  IV-A-4),  they  imply  that  the  material  shocked 
to  states  between  50  and  130  kb  relaxes  along  stresa-particle  velocity 
curves  which  are  nearly  those  which  would  be  expected  if  the  material 
produced  by  the  ahock  were  solid  quartz. 

Above  65  kb  the  Hugoniot  for  Coconino  sandstone  lies  at  smaller  vol¬ 
umes  than  that  of  either  single  or  polycrystalline  quartz.  The  implication 
of  this  behavior  in  relation  to  the  temperature  achieved  during  ahock  com- 
presaion  is  considered  below  (Section  IV-A-6). 

3.  Calculation  or  Shoce  Temperature  in  Sandstones 

The  temperatures  achieved  in  nonporous  rocks  and  minerals  during  shock 
compression  are  rather  moderate.  For  many  nonporous  solids  the  temperature 
achieved  upon  a  shock  compression  and  upon  subsequent  adiabatic  relaxation 
may  be  calculated  using  the  method  described  by  Walsh  and  Christian.1*  To 
date,  temperatures  reached  by  shock  wave  compression  of  porous  rocks  have 
not  been  specifically  discussed.  Temperature  calculations  in  poroua  metals 
are  briefly  discussed  by  Kormer  et  al.a 

Consider  the  shock  compression  of  a  porous  rock  initially  in  a  state, 
given  by  F#,  cr#  ■  0,  to  a  state  V,  c r  (Fig.  12).  From  Eq.  (3)  the  difference 


♦ 


?  **?,**?**  **  'i”1*  e!rr,.r,e,r*  (Fig.  10)  tlM  k«*t  nttniutin  of  tka  pracaraar  aaplitafo 

for  all  thraa  aaaAataaaaliaa  batoaaa  4.1  aaf  $.1  kb  (Tahla  I). 


Pravioaa  tao-Aiaaaaiaaal  ooApa  aa4  plaaa-wara  axpariaaati 
Iha  aarliar  raaalta  aaiag  tha  too-Aiaaaaioaal  aafga  (Raf. 
M)  ara  aoo  kaooa  to  ho  imttrrtel  *ora  IS  kh. 


Sara  atroafly  AiacorAaat  Aaaaitiaa  abova  40  kb. 
I,  Fi|i.  II  aa4  10,  aa4  Tabla  IV,  paftl  IS  awl 
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in  internal  energy  per  unit  maaa  across  the  shock  wave  interface  is 
given  by 


E  -  cr(V0  -  F)/2 


The  difficulty  in  calculating  the  temperature  in  a  shock  transition  lies 
in  determining  how  E  ia  partitioned  into  compressional  energy,  denoted 
by  Ek,  and  thermal  energy,  denoted  by  E 

Thus, 


F, 


E.  +  E. 


(?) 


It  is  convenient  to  discuss  the  energies  of  shock  compression  in 
terms  of  three  regimes  on  the  Hugoniot  curve  'Fig.  12). 


FIG.  12  HUGONIOT,  POROUS  ROCK:  I,  elastic 
regime;  II,  collapse  regime;  QL  energy 
partition  regime.  Shaded  area,  error  in 
energy  from  greater  compression  of 
porous  material 

The  elastic  regime,  or  regime  1,  extends  from  (1'0,0)  to  (1*  ,a  )  where 
the  subscript  e  refers  to  the  Hugoniot  elastic  limit.  For  a  shock  in 
which  the  highest  stress  attained  is  below  c r  ,  the  temperatures  may  be 
estimated  if  adiabatic  compression  ia  assumed.  This  ia  justified  for  the 
low-level  elastic  precursor  observed  in  porous  sandstones  and  in  lime¬ 
stones,  since  the  precursor  shock  velocity  is  close  to  the  measured  adia¬ 
batic  longitudinal  wave  velocity.  The  Hugoniot,  therefore,  may  be  assumed 
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to  lie  close  to  the  adiabat.  The  adiabatic  heating  is  thus  estimated 
from 


-  V) 

Ti t\j  ■  Tt  exP  - Y -  (8> 

where  Tff  is  the  temperature  along  the  Hugoniot,  TQ  is  the  initial  temper¬ 
ature  and  y  is  Gruneisen's  ratio  given  by 


For  several  nonporoua  solids  y/V  is  nearly  constant;14  therefore,  because 
of  the  small  variation  of  ar  and  V  even  for  porous  media  in  regime  I,  the 
use  of  Eq.  (8)  appears  justified.  For  Coconino  sandstone  the  temperature 
increase  achieved  upon  shock  compression  to  the  Hugoniot  elastic  limit, 
from  Eq.  (8),  is  approximately  15°C.  This  temperature  increase  can  be 
neglected  when  it  is  compared  to  the  much  greater  temperature  increase 
resulting  from  stronger  shocks. 

The  second  regime  may  be  defined  as  consisting  of  compressed  states 
lying  between  volumes  V f  and  V  .  In  this  region  the  shock  compression 
may  be  considered  as  primarily  consisting  of  crushing  porous  sandstone  to 
the  volume  of  uncompressed  single-crystal  quartz.  The  inherent  irreversi¬ 
bility  of  such  a  process  suggests  that  nearly  the  entire  energy  of  the 
shock  will  go  into  thermal  energy.  With  this  simplifying  assumption,  the 
recoverable  mechanical  energy  and  contributions  to  E  from  the  quartz-grain 
surface  energy  are  neglected.  The  shock  energy  may  then  be  written  as 

E l„  -  E,  -  cr(V9  -  V)/2  (10) 

For  regime  II  in  the  case  for  silicate  minerals,  the  Debye  temperature  is 
higher  than  the  shock  temperature;  therefore  a  constant  value  of  the  specific 
heat— such  as  6  calorie/mol.  deg.,  the  Dulong-Petit  value — is  inappropriate  to 
calculate  TK  from  E f.  In  quartz,  typical  of  the  common  silicates,  the  specific 
heat  at  constant  volume,  ep,  varies  by  a  factor  of  1.6  between  room  temperature 
and 800°K  (Fig.  13).  It  should  also  be  noted  that,  as  in  the  calculation  of  shock 
temperature  in  porous  metals,13  the  specific  heat  of  the  porous  rock  or  mineral 
is  gi ven  in  terms  of  mass  rather  than  volume  and  is  assumed  to  be  the  same  as 
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FIG.  13  SPECIFIC  HEAT  vs.  TEMPERATURE,  Si02  (a-Quorti) 
(Data  from  Kolloy,  1949) 


that  of  the  solid material ;  thia  being  conaiatent  with  the  •■sumptions  used  in 
obtaining  Eq.  (10).  Preaatire  affects  specific  heat,  within  a  rough  approxima¬ 
tion  discussed  by  Birch,  by  an  amount  proportional  to  the  square  root  of  the 
density  ratio  of  compressed  to  uncompressed  aolid  material.14  In  the  following, 
therefore,  the  specific  heat  is  assumed  to  be  independent  of  stress,  since  the 
compression  of  the  solid  quarts  is  small.  Thus  the  speci  f ic  heat  i s  assumed  to 
be  only  a  function  of  the  temperature  achieved  in  shock  compression. 

At  constant  pressure  cf  is  related  to  c,  by 


c 


r 


(11) 


where  a  is  the  volume  coefficient  of  expansion. 

The  measurements  of  Kelley1*  provide  as  ■  function  of  temperature. 

It  is  then  possible  to  calculate  cy  as  a  function  of  tempersture  (Fig.  13), 


2S 


although  the  difference,  generally,  is  negligible.  Often  c(  or  cf  are 
expressed  as  a  function  of  temperature  in  the  form 


Cv  -  o  +  bT  -  c/r* 


where  a,  6,  and  c  are  experimentally  determined  constants. 


Then 


(12) 


(13) 


For  states  in  regime  II  (Fig.  12)  the  temperature  may  be  implicitly 
expressed  as 


0 


n 


(14) 


where  Eqa.  (12)  and  (10)  were  substituted  into  Eq.  (13).  Eq.  (14)  is  a 
cubic  equation  and  may  be  solved  for  Tt,  the  temperature  achieved  in  any 
shock  compression  to  volumes  in  the  range 

V  <  V  <  V 

»  • 

In  regime  III  (Fig.  12)  a  partition  of  shock  energy  is  assumed.  The 
total  energy  of  a  shock  may  be  written  as 

*lni  "  <*o  ‘  V)cr/2  -  r*C.9(T)dT  +  F.k  (15) 

r» 

The  low  value  of  the  thermal  expansion  coefficient  for  nonmetsllic 
materials  such  as  silicates  and,  in  particular,  quarts,  suggests  thst  the 
0°K  isotherm  may  be  approximated  by  the  isothermal  compression  curve,  such 
as  the  one  determined  statically  by  Bridgman.*  The  0°K  isotherm  msy  be 
expressed  in  terms  of  a  Murnaghan  equation  of  state.  For  this  equstion 
of  state,  the  pressure- volume  relstion  is  written  in  the  form 


where  /?  and  g  are  empirical  conatanta.  The  compreaaional  ahock  energy 
ia  given  by 


E.  -  /  ~  [PdF] 


v 

f 


(17) 


which,  upon  substitution  of  Eq.  (16)  and  integration  becomea 


*(V  /V)*-'  -  1 
*  ■  *,*[ - — -  +  v!\  -  \ 


(18) 


Therefore,  the  temperature  of  a  porous  aandatone  at  final  states  of 
shock  compression  in  regime  III,  may  be  written  implicitly  as 

.  2a  .  sfC'.-n  f iVq/Vf-l-l  1 

•  •  -  >J 


ir’ 

*~r  * 


t) 


Tn  +  2c/b 


(19) 


where  Eq.  (12)  and  Eq.  (17)  are  substituted  into  Eq.  (15).  The  form  of 
Eq.  (19)  is  similar  to  Eq.  (14). 

The  ahock  temperatures  for  Coconino  sandstone  were  cslcnlsted  using 
Eqs.  (14)  snd  (19)  snd  sre  tabulated  in  Table  II.  These  results  snd  the 
temperature  obtained  by  Wackerle  for  quartz  using  the  Welsh  snd  Christisn 
development  are  shown  in  Fig.  14. 

The  temperatures  shown  in  Tsble  II  should  be  considered  as  being 
minimum  temperatures  for  the  following  reasons: 

(1)  It  is  implicitly  sssumed  in  regimes  II  and  III  that 

thermal  equilibrium  behind  the  shock  hes  been  achieved; 
i.t.j  the  shock  energy  density  is  evenly  distributed 
within  the  shocked  msterisl.  In  s  cosrse  sggregste, 
more  thermsl  energy  should  be  concentrsted  at  the  grain 
boundaries  of  the  shocked  meterial  than  in  the  volume 
within  the  center  of  the  quartz  grains,  the  chsrscter' 
istic  relsxstion  time  being  determined  hy  the  grsin 
size  and  thermal  properties  of  the  shocked  material. 


ST 


Table  II 

SHOCK  COMPRESSION  TEMPERATURES,  COCONINO  SANDSTONE 


STRESS 

(kb) 

jr 

(arga/g  *  1  cr®> 

*K 

(arga/g  *  10 "®) 

(arga/g  *  10**) 

T 

(“K) 

0.0 

0.0 

0.0 

0.0 

300 

2.4 

0.0119 

0.0 

0.0119 

305 

0.8 

0.0749 

0.0 

0.0749 

310 

8.3 

0.145 

0.0 

0.145 

322 

12.  S 

0.374 

0.0 

0.374 

348 

17.0 

0.722 

0.0 

0.722 

390 

22.0 

1.20 

0.0 

1.20 

443 

30.0 

2.00 

0.0 

2.00 

524 

33.0 

2. 31 

0.025 

2.28 

556 

39.0 

2.92 

0.15 

2.77 

598 

47.0 

3.76 

0.42 

2.34 

656 

57.5 

4.89 

0.89 

4.00 

713 

73.5 

6.61 

1.61 

5.00 

800 

0.4  0.8  0.8  0.7  as  at  1.0  l.t  1.4  1.8  1.8  t.O 

•K  x  10'*  *»»" 
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(2)  In  general,  the  mineral  in  a  poroua  rock  may  not  remain 
in  the  aame  polymorphic  phaae  aa  the  aolid  material 
ahocked  to  the  aame  preaaure.  In  the  preaent  caae,  for 
Coconino  aandatone,  the  aolid  Hugoniot  of  Fig.  8  appeara 
to  crosa  over  the  curve  (Bridgman* a)  for  the  0°K  iso¬ 
therm  used  in  the  calculation  at  approximately  75  kb. 

An  error  produced  by  too  large  a  value  for  EK  and  there¬ 
fore  too  low  a  temperature  estimate  will  result  from 
using  Eq.  (19);  the  magnitude  of  the  error  in  E K  ia 
equal  to  the  cross-hatched  area  in  Fig.  12. 

(3)  If  the  quarts  in  the  Coconino  sandstone  is  partially 
converted  to  atiahovite  when  shocked  to  above  67  kb, 
as  suggested  in  Section  IV- A- 6  (see  also  Fig.  14),  ad¬ 
ditional  inaccuracies  may  arise  in  the  estimate  of  T H, 
aa  follows: 

(а)  The  °K  isotherm  of  atiahovite  is  of  significantly 
lower  volume  than  that  of  quartz;  thus  an  error 
may  result  by  using  too  high  a  value  for  E K  in 
Eq.  (19).  This,  in  turn,  produces  too  low  a 
value  of  TK. 

(б)  At  the  temperature  of  the  hypothetical  atiahovite 
transition  the  specific  heat  calculated  for 
quartz  has  Tisen  to  that  of  the  classical  Dulong- 
Petit  solid;  an  additional  underestimate  could 
result  if  the  specific  heat  of  atiahovite  were 
actually  lower  than  t  at  of  quartz. 

(c)  If  the  transformation  is  exothermic,  as  implied 
by  the  slope  of  the  quartz-stishovite  phaae  line 
(Fig.  14),  still  another  temperature  underesti¬ 
mate  will  he  made.  For  a  latent  heat  of  2.8  *  18* 
ergs/g  (Section  IV-A-6),  significant  compared  to 
the  thermal  energies  available  above  68  kb  listed 
in  Table  II,  the  phase  transition  itself  could 
provide  an  important  contribution  to  the  shock 
temperature. 

Finally,  it  should  be  noted  that  it  is  not  presently  clear  whether 
the  surface  energy  of  a  porous  rock  will  increase  or  decrease  during 
shock  compression.  Thus,  a  temperature  error  of  unknown  sign  or  magnitude 
resulting  from  surface  effects  may  be  present  in  Table  II. 

4.  Relaxation  Velocity 

The  velocity  of  relaxation  waves  emanating  from  the  lateral  free- 
surfaces  of  shock- compressed  sandstone  specimens  has  been  examined  for 
shock  amplitudes  in  the  60-186  kb  range.  A  study  of  these  unloading 
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velocities  can  provide  information  on  the  properties  of  the  shocked  solid 
which  are  not  obtainable  from  the  Hugoniot. 

For  a  plane  shock  impinging  on  a  specimen  with  a  lateral  boundary 
oriented  at  90°  to  the  shock  front,  the  extent  of  the  region  of  the  inter¬ 
action  of  relaxation  waves  with  the  shock  may  be  determined  by  meaauring 
the  point  at  which  the  planarity  of  the  emerging  wave  is  affected  by  the 
specimen  edge.  The  velocity  of  the  relaxation  wave,  C,  is  related  to  the 
width  of  disturbed  region,  f,  by 

C  -  +  ({/  -  u)2/(  U)2  (20) 

where  D  is  the  thickness  of  the  specimen.  Even  under  ideal  conditions, 
it  is  difficult  to  estimate  W — a  typical  smear  camera  record  is  shown  in 
Fig.  11.  The  problems  which  arise  in  applying  this  method  to  metals  are 
discussed  by  ATtshuler  et  al.u  Due  to  the  inherent  nonplanarity  of  the 
shock  driven  through  a  porous  rock  (Figs.  10  and  11)  the  application  of 
this  technique  to  porous  rocks  especially  at  low  shock  stresses  is  more 
difficult  than  for  metals.  Although  complete  results  were  not  obtained 
for  all  the  experiments  performed,  those  available  are  plotted  in  Fig.  15. 


FIG.  15  LATERAL  RELAXATION  VELOCITY  vs.  MAXIMUM  SHOCK  STRESS 
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The  increase  in  relaxation  velocity  up  to  134  kb  ia  expected.  Above 
that  atreaa  level  the  reaulta  are  aurpriaing.  They  could  indicate  a 
phaae  change  but  thia  apeculation  aeema  at  preaent  inappropriate  in 
view  of  the  inherent  inaccuraciea  in  the  poroua  rock  meaaurementa,  aa 
well  aa  an  inconaiatency  which  ia  noted  below  between  the  lateral  re¬ 
laxation  and  the  free-aurface  velocity  data. 

The  free-aurface  velocity  obtained  in  the  varioua  experimenta  may 
be  compared  to  a  velocity  which  ia  eatimated  from  the  lateral  relaxation 
data.  For  a  aimple  ahock,  the  free-aurface  velocity  is  given  by 
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(21) 


the  first  term  on  the  right  is  the  uaual  particle  velocity  for  shock  com¬ 
pression  from  the  Rankine-Hugoniot  jump  conditions,  and  uf  is  the  relax¬ 
ation  particle  velocity. 

If  the  relaxation  process  is  assumed  to  be  adiabatic,  the  relaxation 
particle  velocity  may  be  written  asn 


Ur  -  rPCdp/p  (22) 

If  C  is  assumed  a  function  of  p  only  (independent  of  temperature) 
the  integral  on  the  right  (Eq.  22)  may  be  evaluated  for  the  data  in 
Fig.  15.  Since  C  has  not  been  obtained  below  60  kb,  the  integral  in 
Eq.  22  is  unknown  for  densities  less  than  2.958  g/cc,  and  the  vertical 
position  of  the  curve  segment  representing  uf  in  Fig.  16  is  undetermined. 
The  point  of  inflection  at  a  density  of  approximately  3.25  g/cc  (Fig.  16) 
corresponds  to  the  apparent  maximum  relaxation  velocity  at  134  kb 
(Fig.  15).  From  comparing  Figs.  15  and  16  it  is  evident  that  the  value 
of  the  integral  in  Eq.  (22)  is  relatively  insensitive  to  the  possible 
large  errors  in  the  determination  of  C. 

When  uf  is  calculated  using  Eq.  (21)  the  slope  of  the  uf  versus  a 
curve  predicted  by  the  free-aurface  velocity  is  in  disagreement  with  that 
obtained  from  the  lateral  relaxation  velocity.  The  factors  producing 
this  apparent  discrepancy  are  not  at  present  known;  however  it  is  possible 
that  aside  from  an  error  in  the  experimental  technique,  which  cannot  be 
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FIG.  16  RELAXATION  PARTICLE  VELOCITY  FROM  LATERAL  AND 
FREE-SURFACE  VELOCITY,  COCONINO  SANDSTONE 


completely  ruled  out,  two  of  the  assumptions  in  the  foregoing  treatment 
may  be  invalid.  These  are: 

(1)  The  shock-compressed  rock  into  which  relaxstion  waves 
propagate  may  be  anisotropic  with  respect  to  the  di¬ 
rections  parallel  and  perpendicular  to  plane  of  the 
shock.  In  consideration  of  the  anisotropy  exhibited 
by  single  crystal  quarts  above  shock  stresses  of  200  kb, 
it  seems  reasonable  to  suppose  that  the  relaxation  path 
may  actually  differ  for  lateral  relaxation  waves  and 
free-surface  reflection  relaxation  waves. 

(2)  The  lateral  relaxation  velocity  may  be  strongly  dependent 
on  shock  temperature  rather  than  being  only  dependent  on 
the  final  compressed  density;  in  this  case  the  application 
of  Eq.  (26)  would  be  invalid  since  relaxation  could  pro¬ 
ceed  along  a  different  path  for  each  temperature. 
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In  view  of  these  and  other  uncertainties,  it  seems  imprudent  to 
speculate  on  this  comparison— at  least  until  more  data  are  available. 


S.  Water  Saturated  Sandstone 

Three  exploratory  experiments  at  high-shock  stress  on  the  effect  of  water 
saturation  on  the  equation  of  state  of  porous  sandstone  (Table  I  and  Fig.  8) 
have  given  points  in  a  fair  agreement  with  the  stress- volume  curve  obtained 
by  Schall1®  using  a  flash  X-ray  technique.  The  rocks  were  saturated  by  means 
of  first  subjecting  them  to  vacuum  and  then,  using  water,  increasing  the 
ambient  pressure  to  atmospheric  pressure.  On  the  basis  of  the  uniformity 
of  shock  velocity,  saturation  was  believed  to  be  essentially  complete. 

At  shock  stresses  sufficiently  high  to  overdrive  the  Hugoniot  elastic 
limit,  no  major  experimental  difficulties  arise.  The  impedance-match 
technique  is  applied  to  obtain  the  equation  of  state  point. 

At  low  stress  levels,  where  a  two-wave  shock  is  likely,  the  use  of 
the  inclined  mirror  technique  (Section  III-A.)  for  measuring  free-surface 
velocity  does  not  seem  appropriate  for  the  following  reasons: 

(1)  To  avoid  observation  of  an  air  shock,  inclined  mirror  measure¬ 
ments  are  usually  performed  in  a  vacuum;  and 

(2)  It  seems  apparent  that  when  a  shock  wave  impinges  on  the  free- 
surface  of  a  saturated  rock,  the  water  exposed  on  the  surface 

and  the  quarts  grains  in  the  sandstone  will  move  with  different  psrticle 
velocities  appropriate  to  the  ahock  impedance  of  each  of  the  components. 

Because  of  these  difficulties,  two  of  the  low-shock  stress  experi¬ 
ments,  one  with  saturated  Coconino  sandstone,  and  the  other  unsaturated,  were 
performed  using  the  quarts  transducer  technique  (Section  III-C).  The 
pressure-time  profiles  of  the  shock  transitions  (Fig.  17)  are  believed  to 
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FIG.  17  QUARTZ  TRANSDUCER  STRESS-TIME  PROFILES,  COCONINO  SANDSTONE 
Nominal  stress  scale  refers  to  rock;  (o)  dry;  (b)  woter  saturated 
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be  quantitatively  reliable  aa  far  aa  the  inherent  timea  in  the  ahocks 
are  concerned.  The  anomaloualy  low  valuea  of  the  precuraor  amplitude 
obaerved  for  dry  sandstone  as  compared  to  thoae  obtained  using  the  in¬ 
clined  mirror  method  (Table  I),  which  are  thought  to  be  the  correct 
values,  indicate  a  discrepancy  presently  unexplained. 

The  elastic  precursor  rise-time  of  0.21  /isec  (Fig.  17),  almost  a 
factor  of  10  greater  than  can  be  attributed  to  tilt  or  nonplanarity  of 
the  shock,  is  a  particularly  interesting  result  for  both  experiments. 

The  quartz  transducer  was  in  the  nonlinear  stress-current  range 
for  the  stress  level  of  the  main  shock  in  both  experiments;  however,  the 
CFO  records  imply  that  the  rise  timea  of  the  main  shocks  are  also  approxi¬ 
mately  0.22  fiatc.  A  rise  time  of  the  main  shock  of  the  order  of  0.2  /zsec 
has  also  been  observed  using  the  inclined  mirror  cutoff  (Figs.  10  and  23), 
in  both  the  sandstone  and  limestone  experiments.  The  precursor  rise- time 
is  difficult  to  observe  in  a  smear  camera  record,  although  its  presence 
may  account  for  some  of  the  nonreproducibility  in  the  measurement  of  the 
precursor  velocity. 

6.  Discussion,  Snock  Compression  op  Sandstone 

The  behavior  of  the  Hugonitt  for  Coconino  sandstone  above  67  kb  is 
particularly  interesting.  Normally  the  Hugoniot  of  the  porous  material 
would  be  expected  to  lie  toward  greater  volume  states  than  that  of  the 
isotherm  and  Hugoniot  of  the  solid  material  because  of  the  heat  generated 
and  the  entropy  gained  in  shock  transition.  In  contrast,  for  shocks 
greater  than  67  kb,  the  volume  achieved  by  the  shock  transition  in 
Coconino  sandstone  is  less  than  that  achieved  at  the  same  pressure  with 
the  solid  material  or  that  observed  in  isothermal  static  compression. 

The  locus  of  pressure-temperatures  states  shown  in  Fig.  14  implies 
that  at  approximately  68  kb  and  770°K  the  stishovite  stability  regime  **is 
achieved  by  shock  compression.  Because  of  what  appears  to  be  essentially 
coincidence  in  shock  stress  between  the  predicted  quartz  (or  coesite?)— 
stishovite  transition  and  an  apparent  anomalous  density  increase  observed 
on  the  Hugoniot,  the  hypothesis  that  stishovite  is  produced  is  quite  sug¬ 
gestive.  The  close  agreement  in  stress  levels  is  undoubtedly  fortuitous. 

A  parallel  hypotheais  has  been  put  forward  by  McQueen  et  al.a  for  the  formation  of 
stishovite  at  144  kb  and  523°K  in  lackerle's  shock  compression  experiment  on  quartz. 

• 

Th*  Sclir  *1  at.  pkaaa  liaa  (Fig.  14)  say  aaffaat  too  kipk  aa  aqaiiibriaa  taaparatara  aaS  praa- 
aara  for  atiakoaita  aiaea  tka  atarkiag  aatarial  *a»  ailica  pal  ratkar  tkaa  eaaaita  ar  qaarta. 
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The  dashed  curve  of  Fig.  8  is  a  possible  configuration  of  the  Hugoniot 
which  is  compatible  with  the  data  points  and  implies  a  polymorphic  tran- 
aition  at  67  kb.  The  polymorphic  transition  indicated  by  the  dashed  curve 
(Fig.  8)  is  not  likely  to  give  rise  to  a  multiwaved  shock  (see  Section  IV-B-2). 

The  stishovite-quartz  energy  difference  at  the  stress-level  at  which 
the  transition  may  have  been  partially  achieved  in  the  present  experiment 
can  be  estimated  from  the  Sclar  et  al.a  results  and  the  pressure- volume  curves 
for  quartz  and  stishovite  assumed  by  McQueen  et  al.  The  internal  energy 
difference  is 

E  -  E  -  T(St  -  S  )  -  a(V  -  V  )  (23) 

where  the  subscripts  s  and  q  refer  to  stishovite  and  quartz,  and  S  ia  entropy. 
Applying  the  Clausiua-Cl apeyron  relation  to  the  Sclar  et  al.  ,  phase  line,  the 
entropy  change, 

S,  -  St  -  -3.6  x  106  ergs/( °K)g 

The  volume  change, 

Vt  -  Vf  -  -0.108  cc/g 

ia  obtained  from  McQueen  et  al.,  and  may  then  be  used  to  estimate  the  energy 
difference  in  Eq.  (23).  For  T  ■  770°K  and  a  •  68  kb, 

Et  -  Ef  -  4.6  x  10*  ergs/g  (24) 

which  may  be  compared  with  the  value 

•  1.5  x  101#  ergs/g 

given  by  McQueen  et  al.  for  the  energy  difference  at  cr  •  0.  In  the  McQueen 
estimate  (a),  the  quartz-stishovite  phase  line  is  assumed  to  pass  through 
the  point  (cr  ■  144  kb,  T  *  523°K),  at  which  the  Wackerle  data  indicates  a 
transition,  and,  (61  the  first  term  in  Eq.  (23)  is  negligible.  The  validity 
of  these  assumptions  and  those  used  in  the  present  calculation  could  be 
critically  considered  when  more  SiOg  phase  equilibrium  data  are  available. 


B.  CARBONATE  ROCKS 

1.  Calcite 

Much  of  this  information  is  taken  from  Semiannual  Report  No.  3, 

October  20,  1962.®  The  calcite  specimens  come  from  naturally  occurring 
“  Iceland  Spar”  found  in  Chihuahua ,  Mexico.  Only  crystals  that  are  colorless, 
free  from  cleavage,  parting,  and  obvious  defects  are  used.  The  specimens, 
usually  about  !^*inch  thick  and  VA  inches  on  a  side,  are  cut  into  four  different 
orientations:  X,  Y,  Z,  and  (0112)-cleavage  plane.  'fhe  designation  of  a  particular 
cut  means  that  the  shock  travels  parallel  to  the  designated  crystallographic  axis. 

A  plot  of  the  experimental  data  in  the  stress- volume  plane  is  shown 
as  Fig.  18;  the  corresponding  stress-particle  velocity  plot  is  shown  as 
Fig.  19.  In  Fig.  20,  hydrostatic  data  reported  by  Bridgman10  and  both 
static  and  dynamic  data  reported  by  Adadurov,  et  al .  }  are  shown. 

The  Hugoniot  elastic  limit  lies  between  18.  5  and  23. 7  kb  for  the  cuts  in¬ 
vestigated —  a  conclusion  based  on  more  than  twenty  measurements.  This  tran¬ 
sition  is  termed  the  Hugoniot  elastic  limit  because  it  is  the  first  observed 
and  the  velocities  with  which  the  wave  front  travels  in  various  crystallographic 
directions  correspond  within  3%  to  elastic  wave  velocities  computed  from 
elastic  constants.  Thus  it  would  seem  that  the  value  of  20  kb  is  the  maximum 
normal  stress  that  calcite  can  support  without  exceeding  the  maximum  values  of 
lateral  shear  stresses  that  can  exist  in  calcite,  t .  e. ,  a  true  Hugoniot  elastic 
limit.  But  this  may  not  be  so.  In  Vermont  and  in  Yule  marble  a  Hugoniot  elastic 
limit  of  about  12  kb  is  measured  and  a  transition  at  22  kb.  Hence  maximum 
elastic  amplitudes  in  calcite  may  be  controlled  by  a  polymorphic  transition 
at  22  kb  rather  than  being  limited  by  the  maximum  lateral  shear  stresses 
that  calcite  would  otherwise  support. 

Good  evidence  exists  in  al  1  orientations  for  a  phase  transition  at  30 
to  32  kb;  about  eight  measurements  are  available.  The  transi tion  at  95  kb 
is  also  wel  1  documented;  about  four  measurements  are  available  for  all 
orientations  in  calcite.  The  transition  at  45  kb  is  least  documented;  only  two 
reliable  values  are  avai lable  in  Z-cut  calcite,  with  some  evidence  that  there 
is  a  transition  at  about  45  kb  in  three  other  calcite  experiments. 

Dremin  and  Adadurov21  record  a  transition  in  marble  at  150  kb.  Their 
technique  measures  only  an  average  free-surface  velocity;  i.e.,  above 
150  kb  the  free-surface  velocity  of  a  single  wave  front  was  measured, 
and  below  150  kb  the  average  free-surface  velocity  of  a  multiple  wave 
front  was  measured.  We  did  not  observe  a  transition  at  150  kb  but  did 
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FIG.  19  STRESS-PARTICLE  VELOCITY  PLOT  FOR  CALCITE 
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FIG.  20  STRESS*p0/p  PLOT  FOR  CALCITE 
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observe  one  at  95  kb.  The  150  kb  transition  recorded  by  Adadurov  et  al.  ,* 
is  probably  the  boundary  between  a  double  and  single  shock  front  produced 
by  the  95  kb  transition.  If  Rayleigh  lines  are  constructed  along  the 
calcite  Hugoniot,  one  can  see  that  a  single  wave  is  stable  for  pressure 
states  above  150  kb  and  a  multiple  front  is  stable  below  this  region. 

The  complicated  shock  front  below  150  kb  was  not  observed  by  Adadurov 
et  al.,  because  the  gap  technique  does  not  resolve  such  multiple  wave 
fronts. 

The  differences  in  the  data  among  the  various  crystal  cuts  at  250  kb 
are  disturbing  (a  similar  behavior  is  observed  for  the  various  cuts  of 
quartz  at  these  pressures).  They  imply  that  anisotropy  is  still  important. 
It  would  be  desirable  to  repeat  these  measurements  since  there  exists  only 
one  experiment  for  each  cut  at  these  pressures. 

2.  Marbles  and  Limestones 

Yule  and  Vermont  marble,  Solenhofen  and  Spergen  limestones  were  used 
in  the  experiments  to  measure  Hugoniot  data  for  polycrystalline  carbonate 
rocks.  A  petrographic  description  of  these  rocks  follows: 

a.  SOLENHOFEN  LIMESTONE 

Solenhofen  limestone  is  a  massive,  tan,  lithographic  limestone 
composed  of  96%  calcite. 

The  grain  size  of  the  calcite  ranges  from  0.005  to  0.015  mm  in 
diameter  and  has  a  median  diameter  of  0.010  mm.  There  are  occasional 
grains  as  large  aa  0.03  mm  and  an  occasional  quartz  grain.  A  few  percent 
of  interstitial  clay  is  present  and  the  true  porosity  is  a  few  tenths  of 
a  percent.  The  density  is  2.581  g/cc. 

b.  SPERCEN  LIMESTONE 

The  Spergen  limestone  is  a  gray,  well  cemented,  microcoquina 
composed  of  about  60%  fossil  fragments,  20%  bioclastic  debris  including 
carbonate  detritus,  20%  calcite,  and  a  percent  or  two  of  oolites. 

The  size  of  the  sub-rounded  fragments  ranges  from  0.1  to  0.7  mm 
with  the  average  size  being  about  0.3  to  0.4  mm  in  diameter.  The  average 
grain  size  of  the  calcite  matrix  is  0.2  to  0. 3  mm  in  diameter  with  an 
occasional  grain  as  large  as  1.0  mm  in  diameter. 
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The  calcite  matrix  tends  to  have  a  “ sub-poiki loblastic ”  texture. 
Some  of  the  bioclastic  fragments  which  are  included  in  the  grains  of 
cementing  calcite  have  the  same  optic  orientation  as  the  neighboring 
calcite  grains. 

The  porosity  is  about  12%.  The  density  ranges  from  2.360  to 
2.398  g/cc. 


C.  VERMONT  MARBLE 

Vermont  marble  is  massive,  white,  fine-grained  and  composed 
of  95%  calcite  grains. 

The  calcite  ranges  from  0.05  to  0.2  mm  in  diameter  and  has  a 
median  diameter  of  0.12  mm.  Five  percent  secondary  quartz  occurs  in 
occasional  bands  that  are  1  mm  thick.  The  quartz  grains  are  0.05  mm  in 
diameter  and  are  scattered  among  the  calcite  grains  in  the  1  mm  band. 

The  density  is  2.687  cc/gm. 

d.  YULE  MARBLE 

Yule  marble  is  a  massive,  white,  homogeneous  marble  that  is 
composed  of  100%  calcite. 

The  average  grain  size  of  the  marble  is  between  0.2  and  0.4  mm 
with  an  occasional  grain  as  large  as  1.5  mm  in  diameter  (apparent  grain 
size  in  the  thin  section  depends  upon  orientation  as  the  grains  are  platy). 
The  density  is  2.697  cc/gm. 

Table  III  compiles  all  plsne-wave  experimental  data  on  carbonate 
rocks,  including  a  reinterpretation  of  the  plane-wave  experiments  listed 
in  Semiannual  Report  No.  2.  Figure  18  is  a  plot  of  the  llugoniot  curves 
for  the  carbonates. 

The  Hugoniot  for  marble  follows  the  llugoniot  for  calcite.  The 
Hugoniot  elastic  limit  (10-15  kb)  and  3  transitions  (22,  40,  and  95  kb) 
are  observed. 

For  Vermont  marble  the  Hugoniot  elastic  limit  is  about  10  to 
15  kb  and  decays  as  the  shock  travels.  Figure  21  shows  Hugoniot  elastic 
limits  plotted  as  a  function  of  specimen  thickness.  Data  for  '«>inch  and 
%-inch  specimens  of  Vermont  marble  are  measured  from  a  single  experiment 
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FIG.  21  DECAY  OF  HUGONIOT  ELASTIC  LIMITS 
Final  Stata  Pressure*  art  Indicated 
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and  are  conaidered  reliable  valuea.  Data  for  Z-cut  calcite  are  taken 
from  several  experiments.  In  general,  both  sets  of  data  show  that  the 
Hugoniot  elastic  limit  decays  at  about  the  same  rate. 

The  experiment  using  a  J^-inch  and  a  %-inch  specimen  of  Vermont 
marble  also  provides  information  about  the  behavior  of  transitions  above 
the  Hugoniot  elastic  limit.  First  consider  what  shock  theory  would  predict. 

A  shock  front  may  become  unstable  if  a  phase  transition  is 
present  in  the  material.  The  phase  transition  may  be  either  a  first- 
order  transition  having  a  discontinuous  change  in  volume,  or  a  second- 
order  transition  having  a  discontinuous  change  in  heat  capacity,  isothermal 
compressibility,  or  coefficient  of  thermal  expansion.22  If  a  shock  front 
becomes  unstable  it  will  break  up  into  two  or  more  shock  fronts,  each 
stable  for  a  certain  range  of  pressures. 

The  necessary  and  aufficient  condition  for  stability  of  a  single 
shock  transition  between  two  arbitrary  states,  a  and  /3  ahould  not  interaect 
the  Hugoniot  at  any  intermediate  point.  If  the  Hugoniot  ia  such  that  the 
Rayleigh  line  does  intersect  it  at  intermediate  statea,  then  the  Rayleigh 
line  must  consist  of  segments,  each  segment  corresponding  to  a  wavefront 
(Fig.  22).  For  a  material  with  several  phase  transitions,  therefore,  the 
wave  structure  can  be  somewhat  complicated  and  can  exhibit  one  or  more 
wavefronta  depending  on  the  final  pressure. 

The  intersection  of  the  Rayleigh  line  with  the  Hugoniot  indicatea 
the  pressure  of  the  transition  or  final  state,  while  the  slope  of  the 
Rayleigh  line  is  a  meaaure  of  the  velocity  for  a  particular  shock  trsnsi- 
tion.  In  fact,  it  ia  easily  shown  that  shock  velocity  is  proportional  to 
the  square  root  of  the  negative  slope  of  the  Rayleigh  line. 

Figure  22  illustrates  a  Hugoniot  for  a  material  with  two  phase 
transitions  in  addition  to  a  yield  point.  For  final  shock  pressures 
greater  than  /,  a  single  shock  is  stable.  For  final  pressures  between  I 
and  E,  or  between  B  and  D,  two  wavea  will  be  present;  the  first  is  the 
elastic  precursor  wave  whose  amplitude  is  defined  ss  the  Hugoniot  elastic 
limit,  and  the  second  is  the  final  state.  For  final  pressures  between  D 
and  E,  three  waves  will  be  present;  the  first  is  the  elastic  transition, 
the  second  is  a  phase  transition,  and  the  third  is  the  final  state.  Note 
that  the  transition  at  C  does  not  give  rise  to  a  separate  wavefront. 

Note,  too,  that  the  region  near  Rayleigh  line  BDE  would  indicate  the 
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transition  at  D  but  it  will  take  tine  until  the  fronts  have  traveled  and 
separated  sufficiently  for  the  recording  equipment  to  resolve  the  fronts. 

For  regions  in  which  there  are  finite  rise  times  to  fronts,  or 
for  regions  near  the  transition  at  D  where  the  two  fronts  may  not  be  coni' 
pletely  resolved,  the  inclined-mirror  trace  appears  as  in  Fig.  23.  An 
inclined-mirror  trace  for  a  steep  front  appears  as  in  Fig.  24. 

For  an  inclined-mirror  trace  from  a  steep  fronted  shock,  the 
arrival  times  and  changes  in  trace  angles  are  abrupt,  but  this  is  not  so 
for  shocks  that  have  a  finite  rise  time;  the  inclined-mirror  traces 

46 


appear  curved.  The  Rankine-Hugoniot  jump  conditiona  are  atill  valid  in 
regions  between  initial  and  final  states  of  a  shock  with  a  finite  rise 
time  and  are  used  piecewise  to  obtain  states  between  the  final  and 
initial  states.  Figure  23  shows  how  the  curved  inclined-mirror  trace 
was  approximated  by  a  series  of  shock  steps.  The  Rankine-Hugoniot  jump 
conditions  was  applied  to  each  step.  In  this  case,  the  states  between 
the  initial  (no.  1)  and  final  (no.  2)  in  Fig.  23  are  non-ateady  atatea 
and  describe  the  path  from  the  first  state  to  the  final  state.  Two  of 
these  paths  are  plotted  in  Fig.  18. 

We  stated  that  additional  information  is  available  for  transi¬ 
tions  above  the  Hugoniot  elastic  limit.  The  inclined-mirror  trace  for 
Vermont  marble  is  shown  in  Fig.  25.  The  trace  from  the  K-inch  specimen 
shows  three  distinct  breaks  indicating  the  front  is  composed  of  three 
shocks.  The  trace  from  the  %-inch  specimen  shows  breaks  which  are  not 
as  distinct  and  are  curved — indicating  that  the  front  is  composed  of 
fewer  shocks  of  longer  rise  times.  Figure  26  is  a  diagram  of  these  two 
fronts.  The  states  in  the  'A-  and  Ji- inch- specimen  are  close  to  the 
Rayleigh  line  path  BDE  of  Fig.  22  thua  it  takes  some  time  for  the  shocks 
to  separate  sufficiently  to  show  the  transitions.  Notice,  that  the 
Hugoniot  elastic  limit  is  decaying;  however,  the  other  two  transitions 
do  not  appear  to  decrease  appreciably— the  observed  variations  are  not 
experimentally  significant. 

Porous  carbonates  compact  irreversibly  and  follow  the  marble 
calcite  Hugoniot  once  they  have  compacted.  Except  for  the  elastic 
transition  at  5  kb,  other  transitions  are  not  observed  because  the  thermo¬ 
dynamic  path  followed  by  the  locking  carbonate  does  not  permit  such  reso¬ 
lution.  This  is  a  similar  case  to  that  part  of  Fig.  22  described  by  the 
Rayleigh  line  BH  near  the  transition  at  G.  The  two  dashed  lines  show 
such  a  particular  path  through  which  the  porous  carbonates  compress  to 
the  final  state.  Note  that  the  Solenhofen  limestone  behaves  as  if  it 
were  a  pure  carbonate  rock  with  a  4%  porosity.  This  is  the  first  indica¬ 
tion  that,  in  aome  caaes,  minor  impurities  in  the  rocks  may  have  con¬ 
siderable  effect  on  propagation  of  shock  waves. 

In  summary,  then,  we  may  describe  the  Hugoniots  for  carbonate 
rocks  as  follows: 

The  nonporous  polycrystalline  carbonates  follow  the  Hugoniot 
of  calcite.  Transitions  at  22,  42,  and  95  kb  are  observed.  The  transition 
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FIG.  25  INCLINED  MIRROR  TRACE  FOR 
VERMONT  MARBLE 
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FIG.  26  IDEALIZED  SHOCK  FRONT  FOR  MARBLE 
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at  30  kb  was  not  obaerved,  but  a  apecial  experiment  would  have  to  be 
designed  aince  it  ia  the  moat  difficult  to  detect.  The  tranaition  at 
22  kb  in  marble  ia  puzzling  aince  thia  ia  the  preaaure  of  the  calcite 
Hugoniot  elaatic  limit.  The  22  kb  tranaition  waa  alao  obaerved  in  the 
two-dimensional  wedge  data.1  Hugoniot  elaatic  limits  are  about  10  to 
IS  kb  and  decay. 

The  porous  carbonates  have  Hugoniot  elastic  limits  of  about 
5  kb  and  compress  to  states  that  lie  along  the  calcite  Hugoniot. 
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C.  PLAGIOCLASE  ROCKS  AND  BASALT 


1.  Plagioclase  Rocks 

Most  experiments  hsve  used  anorthosite.  This  is  an  igneous  rock 
composed  largely  of  the  andesine  member  of  the  plagioclase  series;  our 
specimens  were  collected  in  the  western  San  Gabriel  mountains  in  southern 
California.  The  specimens  were  light  bluish-gray.  About  98%  of  the  rock 
is  andesine  (An  43),  and  the  remainder  consists  of  apatite,  zircon,  chlo¬ 
rite,  hornblende,  biotite,  and  iron  ore  minerals  scattered  throughout. 

The  andesine  grains  are  anhedral  to  aubhedral  and  range  from  about  8  mm 
to  1  mm  in  diameter;  the  average  grain  size  is  about  4  mm.  The  secondary 
minerals  are  usually  imbedded  in  the  andesine  and  range  from  0.01  mm  to 
0.03  mm  in  size.  Density  is  2.662  g/cc. 

Specimens  of  single  crystal  albite  and  labradorite  were  also  used 
in  one  experiment. 

Table  IV  presents  all  data  on  plagioclase.  Ilugoniota  are  plotted 
in  Fig.  27,  stress-particle  velocity  curves  in  Fig.  28.  Figure  29  dis¬ 
plays  all  known  data  for  rocks  which  contain  large  amounts  of  plagioclase. 

The  main  features  of  shock  propagation  in  plagioclase  are  best  demon¬ 
strated  by  the  Hugoniot  for  anorthosite,  and  the  Hugoniot  for  single 
crystal  albite  and  labradorite,  Fig.  27.  We  see  that  these  rocks  have  a 
Hugoniot  elastic  limit  at  about  40  kb  and, as  yet,  have  no  detectable  tran¬ 
sitions  beyond  this  at  least  up  to  pressures  of  250  kb.  Composition  and 
elastic  wave  decay  govern  the  specific  value  of  the  Hugoniot  elastic  limit 
observed  in  a  particular  experiment. 

In  anorthosite  evidence  for  elastic  wave  decay  (stress  relaxation)  is 
indirect;  the  elastic  wave  amplitude  was  observed  to  increase  as  the  ampli¬ 
tude  of  the  second  wave  increased  for  the  same  specimen  thickness.  This 
effect  was  also  observed  in  quartz  by  Kowles*  and  Wackerle.*  One  might 
speculate  that  if  sufficient  single  crystal  data  for  feldspar  were  available, 
the  data  would  show  elastic  wave  decay  to  be  as  important  in  feldspar  as  in 
quarts  —  a  reasonable  speculation  since  crystal  structure  for  the  two 
minerals  is  similar  in  many  respects. 

We  stated  that  the  Hugoniot  elastic  limit  also  depended  upon  the 
composition  of  the  feldspar.  Evidence  for  this  is  meager:  one  experiment. 
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FIG.  27  HUGONIOT  FOR  PLAGIOCLASE  AND  BASALT 
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FIG.  28  STRESS-PARTICLE  VELOCITY  CURVE  FOR  PLAGIOCLASE 
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However,  the  two  feldspars,  lsbrsdorite  and  albite,  have  essentially  the 
same  thickness  and  crystal  orientation  and  were  shocked  by  the  same  input 
pressure.  We  recorded  a  difference  of  about  20  kb  for  their  Hugoniot 
elastic  limit.  One  would  not  expect  elastic  wave  decay  to  be  markedly 
different  in  comparing  the  adjacent  specimens.  Thus,  it  seems  reasonable 
to  expect  Hugoniot  elastic  limits  to  increase  from  more  sodic  members 
to  more  caicic  plagioclase  members. 

The  feldspars  are  a  complicated  group  of  minerals  but  the  above 
experiments  have  given  us  some  insight  into  their  behavior  under  shock 
loading.  Much  work  should  be  done  to  obtain  Hugoniot  data  for  three 
crystallographic  orientations  for  at  least  two  end  members  of  the 
plagioclase  series  and  orthoclase.  Additional  work  is  needed  to  obtain 
values  of  stress- relaxation  for  plagioclase  of  a  specific  composition. 

2.  Basalt 

Vacaville  basalt  is  massive,  dark  gray  basalt  that  can  be  classified 
as  2312E  in  the  system  of  Johannsen.  The  mineral  compoaition  is  53% 
plagioclase  (An$,-An4$),  31%  augite,  9%  magnetite-ilmenite,  5%  celadonite, 
and  2%  apatite.  The  basalt  is  composed  esaemtially  of  equidimensional 
mineral  grains  with  interstitial  celadonite.  The  feldspars  and  augite 
crystallized  penecontemporaneoualy ,  giving  the  rock  an  intergranular 
texture.  Grain  sizes  are  between  0.2  mm  and  0.02  mm.  The  presence  of 
some  vesicules  gives  the  rock  a  porosity  of  about  3%,  Density  is 
2.82  g/cc. 

Table  V  lists  Hugoniot  data.  Figures  27  and  30  are  plots  of 
Hugoniot  data  and  stress-particle  velocity  data  for  basalt.  Figure  29 
shows  all  available  data  for  feldspar  and  basalt  type  rocks. 

The  Hugoniot  shows  cwo  features  of  interest.  The  first  is  the 
Hugoniot  elastic  limit  at  about  50  kb  with  no  other  transitions  evident 
to  pressures  as  high  as  190  kb.  The  second  is  that  velocities  of  the 
main  shock  are  close  to  elastic  velocities.  This  makes  interpretation 
of  a  multiple  shock  front  difficult  because  of  the  time  required  for 
the  fronts  to  separate.  We  do  not  have  enough  data  to  tell  if  the  elastic 
wave  in  basalt  decays  with  distance  of  travel. 

The  transition  at  about  150  kb  reported  by  Hughes  and  McQueen3  was 
not  found  in  the  anorthosite  or  the  basalt.  Moreover,  a  transition  with 
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FIG.  30  PARTICLE  VELOCITY  CURVE  FOR  BASALT 


magnitude  of  the  large  suggested  volume  chsnge  for  gsbbro  seems  unlikely 
in  the  feldspar  rocks  so  far  studied.  The  properties  of  the  gsbbro, 
which  is  composed  of  more  than  50%  plagioclase,  are  probably  similar  to 
those  of  anorthosite.  A  modified  explanation  is  therefore  proposed  for 
the  data  presented  by  Hughes  and  McQueen.  The  extension  of  the  Kayleigh 
line  above  the  Hugoniot  elastic  limit  intersects  the  llugoniot  curve  at 
approximately  150  and  230  kb  for  anorthosite  and  basalt  (Fig.  27).  It 
seems  reasonable  to  conclude  that  rocks  with  a  high  feldspar  content 
would  have  a  stable  single  shock  front  above  a  shock  stress  on  the  order 
of  150  to  230  kb,  and  a  double  shock  front  below  this  magnitude.  Hughes 
and  McQueen  used  a  technique  that  recorded  the  shock  velocity  of  the 
first  front  of  a  multiple  shock  until  they  exceeded  300  kb  (Fig.  31). 
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FIG.  31  SHOCK-PARTICLE  VELOCITY  CURVE  FOR  FELDSPAR  BEARING  ROCKS 


Above  about  300  kb  a  aingle  ahock  again  becomes  the  stable  front  and 
the  ahock  velocity  of  a  single  front  is  measured.  The  effect  recorded 
by  Hughes  and  McQueen  is,  then,  the  boundary  between  regions  of  s  stable 
single  shock  and  a  stable  double  shock.  Their  method  is  based  on  an 
impedance  match  solution  which,  as  they  point  out,  is  invalid  for  a  two- 
wave  structure.  If  a  two-wave  structure  is  present,  it  cannot  be  ob¬ 
served  directly  but  its  presence  can  be  inferred  from  the  change  in  slope 
of  the  shock  velocity-particle  velocity  curve. 

Using  their  system,  the  transition  point  where  the  second  wave 
begins  cannot  be  precisely  determined.  They  suggest  that  the  transition 
is  at  ISO  kb.  However,  the  first  wave  through  gabbro  may  be  an  elastic 
precursor  of  about  SO  kb  rather  than  the  150  kb  state  suggested.  Although 
the  material  of  the  present  investigation  is  not  identical  it  appears 
likely  that  Hughes  and  McQueen  recorded  the  elastic  wave  until  it  wss 
overdriven  at  300  kb,  and  not  a  polymorphic  transition. 


V  STATIC  MEASUREMENTS 


A.  EXPERIMENTS 

(Portions  of  the  following  pages  have  been  taken  from  another  report 
prepared  at  these  Laboratories  for  Vela-Uniform.9) 

The  goal  of  the  static  experiments  is  to  find  some  parameter  which 
can  be  related  to  the  Hugoniot  elastic  limit.  Some  of  the  rocks  studied 
dynamically  have  been  subjected  to  uniaxial  compression,  the  same  strain 
geometry  that  exists  in  the  dynamic  experiments.  Length,  compressions! 
velocity,  and  shear  velocity  have  been  measured  as  a  function  of  presstre. 
Present  results  indicate  that  strength  under  conditions  of  plane-strain 
may  be  directly  related  to  the  Hugoniot  elastic  limit  for  porous  quarts 
sandstones. 

The  details  of  the  theory,  equipment,  and  techniques  are  described 
elsewhere.9  Briefly,  a  cylindrical  rock  sample  is  placed  in  a  steel  pres¬ 
sure  vessel  (Fig.  32).  Pressure  is  exerted  in  the  axial  direction  by 
means  of  two  steel  pistons  which  have  ultrasonic  ceramic  transducers 
bonded  to  the  ends  with  phenylsalicylate.  Lateral  pressure  is  applied 
through  a  hole  in  the  chamber  wall  and  is  controlled  independently  of 
the  axial  pressure  (Fig.  33). 

The  specimen  is  fitted  into  a  rubber  sleeve  with  the  ends  of  the 
sleeve  wired  firmly  to  ends  of  the  steel  pistons  to  prevent  oil  from  pene¬ 
trating  the  sample.  This  technique  is  quick  and  is  usually  quite  suc¬ 
cessful. 

For  plane-strain,  the  relationship 

rtJ  -  [1  -  2(/0/o)*]rl, 

must  be  satisfied,  where  t8J  is  the  lateral  pressure,  is  the  longi¬ 
tudinal  pressure,  /3  is  the  shear  velocity,  and  a  is  the  compressional 
velocity.  When  this  relationship  is  determined,  the  specimen  can  be 
compressed  under  conditions  of  plane-strain.  In  practice  the  shear 
velocity  is  quite  difficult  to  determine  accurately,  especially  for 
porous  sandstones.  At  low  pressures  both  velocities  are  extremely 


FIG.  32  EXPERIMENTAL  PRESSURE  VESSEL  FOR  ULTRASONIC  MEASUREMENTS 
WITH  TWO-DIMENSIONAL  STRESS  SYSTEMS 
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FIG.  33  SCHEMATIC  OF  SYSTEM 


sensitive  to  saall  changes  in  pressure.  The  plane-strain  condition  is 
best  deterained  froa  Tallies  taken  at  pressures  of  about  0.3  kb  or  aore. 
Care  aust  be  taken  to  aessure  both  coapressionel  and  shear  eelocities 
at  the  ssae  pressure  since  saall  differences  in  Telocity  can  aake  s  large 
difference  in  the  plane-strain  line. 

A  1  ixsec,  30-Tolt  video  pulse  is  applied  to  either  of  two  1  Mc/sec 
ceraaic  transducers.  The  trsnsaitted  signal  froa  the  receiring  trans¬ 
ducer  is  observed  on  the  oscilloscope.  Zero-pressure  travel  tiae  aessure- 
aents  are  obtained  using  the  calibrated  Tektronix  Type  535A  oscilloscope 
delay  lines.  As  pressure  is  applied  to  the  speciaen,  variations  of 

ft 


0.005  i±»ec  in  travel  time  can  be  measured  by  displacement  of  the  signal 
on  the  oacilloacope  face  uaing  a  sweep  apeed  of  0.1  ptaec/cm. 

Longitudinally  polarised  ceramic  tranaducera  are  uaed  to  generate 
and  detect  longitudinal  waves.  Tranaverae  velocitiea  are  meaaured  in 
two  ways:  (1)  A  direct  measurement  can  be  made  uaing  square  ceramic 
transducers,  polarized  parallel  to  their  silvered  faces.  These  trans¬ 
ducers  will  produce  a  strong  transverse  pulse  which  is  preceded  by  a 
weak  longitudinal  forerunning  wave.  (2)  An  indirect  measurement  depends 
on  conversion  of  longitudinal  waves  to  transverse  waves  at  the  specimen 
lateral  boundary:  The  transverse  waves  ( 5 )  propagate  acroas  the  specimen 
diameter  and  are  again  reflected  as  longitudinal  waves  (P) .  When  a  pulse 
is  sent  down  the  specimen  by  one  transducer,  the  received  signal  at  the 
other  transducer  consists  of  the  direct  P  arrival  followed  by  a  train  of 
pulses  corresponding  to  paths  of  the  type  P  -  nS  -  P,  where  n  is  an 
integer.  Hie  delay  time  of  these  pulses  is  given  by 

At,  - 

where  D  is  the  cylindrical  specimen  diameter.  The  indirect  measurement 
has  the  advantage  that  both  a  and  fi  can  be  meaaured  at  the  same  time  on 
the  same  sample,  but  suffers  from  the  disadvantage  noted  by  Hughes  and 
Maurette14  that  it  sometimes  gives  spurious  results  on  coarse-grained  rocks. 
In  the  present  work,  this  method  was  not  successful  using  porous  sandstones 
and  had  limited  success  with  other  rocks.  Another  difficulty  with  the  in¬ 
direct  measurement  is  that  it  is  based  on  the  assumption  that  the  medium 
is  homogeneously  strained,  as  by  hydrostatic  pressure.  In  our  experiment, 
the  stress  has  a  tensorial  character  and  hence  mode  conversion  at  the 
boundaries  becomes  more  complicated. 

B.  RESULTS 

1.  St.  Pbtbn  Sandstone 

Moat  of  the  study  has  been  concentrated  on  St.  Peter  sandstone  in 
an  attempt  to  improve  techniques  and  to  obtain  reliable  dsta.  The  most 
interesting  aspect  proves  to  be  the  strength  of  the  rock  under  conditions 
of  plane-strain  or  near  plane-strain.  In  seven  runs,  the  sandstone  has 
failed  at  6.5  ±  0.3  kb.  The  conditions  at  failure  are  somewhat  offset 
from  plane-strain  because  the  maximum  longitudinal  pressure  obtainable 
along  the  St.  Peter  plane- strain  line  is  about  5.7  kb  with  our  equipment. 
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Our  limiting  factor  is  lateral  streas,  which  cannot  be  raised  much  above 
30,000  psi.  The  ratio  of  ram  pressure  to  lateral  pressure  at  failure  was 
slightly  greater  than  that  calculated  for  plane-strain.*  However,  failure 
is  obtained  for  several  different  values  of  ram  to  lateral  pressure  ratio, 
all  larger  than,  but  close  to,  the  ratio  for  plane-strain.  The  strength 
does  not  vary  significantly  with  the  ratio  over  the  small  range  studied. 

The  conclusion  is  reached  that  the  strength  of  St.  Peter  sandstone  under 
conditions  of  plane-strain  is  close  to  6.5  kb. 

Only  two  values  of  the  Hugoniot  elastic  limit  have  been  measured  in 
shock  experiments  10.9  kb  and  4.5  kb.  The  lower  value  is  thought  to  be 
more  reliable.  The  spread  is  large,  but  the  agreement  with  the  strength 
under  static  conditions  is  encouraging. 

Two  examples  of  specimens  after  failure  are  shown  in  Figs.  34  and 
35.  The  2-inch  samplea  have  a  tendency  to  fail  along  a  45°  plane,  as 
shown.  The  1-inch  samples  tend  to  fail  in  a  conical  manner  with  one  end 
serving  as  the  base  of  a  cone,  again  at  about  45°  to  the  axis.  Presumably, 
the  shorter  length  does  not  allow  a  plane  to  form.  In  the  pictures  the 
rubber  sleeve,  of  double  thickness,  has  been  slit  lengthwise  to  expose  the 
samples  . 

Measurements  of  the  change  in  travel  time  of  the  compressional  wave 
with  increasing  pressure  are  shown  in  Fig.  36.  The  values  change  in  a 
characteristic  and  repeatable  manner,  but  the  magnitude  of  the  change, 
especially  over  the  initial  kilobar  or  so,  fluctuates  from  sample  to 
sample.  The  method  of  analysis  depends  on  correcting  for  times  through 
the  transducers  and  pistons  by  obtaining  the  difference  in  measurements 
on  specimens  of  different  lengths,  usually  2  inches  and  1  inch.  Obviously, 
specimen  to  specimen  variations  are  troublesome  in  this  sort  of  analysis. 

In  the  present  case  the  average  values  for  three  runs  on  1-inch  samples 
have  been  subtracted  from  the  average  of  three  runs  on  2-inch  samples. 

The  dial  readings  (Fig.  37),  which  measure  change  in  length  of  the 
samples  and  pistons,  are  reduced  in  a  similar  manner  to  eliminate  the 
effects  of  changes  in  the  piston.  The  dial  readings  are  subject  to 
errors  in  friction  and  in  alignment  of  the  piston  with  the  surface  of  the 
tie  rod  frame.  In  nearly  all  cases,  the  readings  have  changed  with  pres¬ 
sure  in  a  linear  fashion  except  near  the  beginning  of  the  run  where  they 
deviate  markedly,  apparently  through  alignment  effects  and  friction. 

Aa  apparaat  aacaaa  ia  raw  praaaara  a«y  ba  partially  offaat  by  tba  frictioa  ia  tba  "0*  riapa  abich  aaal 
tba  ataal  piatoaa  (Fip.  32).  Ibia  frictioa  ia  tboopbt  to  ba  avail  aai  aill  taaS  to  ra4aca  tba  actwl 
atraaa,  applio4  to  tba  apaaiaaa. 

64 


FIG.  34  TWO-INCH  SPECIMEN  OF  ST.  PETER  SANDSTONE 
AFTER  FAILURE 


FIG.  35  ONE-INCH  SPECIMEN  OF  THE  ST.  PETER 
SANDSTONE  AFTER  FAILURE 


FIG.  36  ST.  PETER  SANDSTONE  -  RAW  ULTRASONIC  DATA 
FOR  COMPRESSIONAL  WAVE  ARRIVALS 


Because  of  errors  at  the  beginning,  slopes  over  the  main  portion  of  the 
plots  are  averaged  and  the  average  slope  is  extended  back  to  aero  pres* 
sure.  The  Measured  decrease  in  sample  length  shown  in  Fig.  6  is  the  dif¬ 
ference  between  the  average  for  three  runs  with  l*inch  specimens  and  the 
average  for  three  runs  with  2- inch  specimens  obtained  in  this  manner. 
Another  run  with  a  2-inch  sample  has  been  disregarded  because  it  differs 
substantially  from  the  other  three.  The  cause  of  the  discrepancy  is  not 
certain,  although  some  difference  evidently  arises  because  the  dial  gauge 
had  been  aligned  somewhat  differently.  For  future  work  the  system  of 
measuring  sample  lengths  should  be  improved. 

The  length  decrease  of  the  sample  may  be  calculated  independently 
of  the  dial  measurements,  using  the  relationship 
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FIG.  37  ST.  PETER  SANDSTONE  -  CALCULATED  AND  MEASURED 
LENGTH  DECREASE 


where  5  is  the  specimen  length,  S„  is  the  originsl  length,  p#  is  the 
original  density,  tf  is  the  travel-time  through  the  specimen  for  the  com- 
pressional  wave,  and  ^  is  the  longitudinal  stress.  This  equation  is 
derived  elsewhere**  and  is  only  valid  for  conditions  of  plane-strain. 

The  calculated  length  decrease  (Fig.  37)  is  in  poor  agreement  with 
the  length  decrease  as  given  by  the  dial  measurements.  As  noted  above, 
the  dial  measurements  are  subject  to  large  errors  in  this  analysis. 

The  measurements  of  length  and  travel-time  permit  the  compressional 
velocity  as  a  function  of  pressure  to  be  plotted.  The  results  are  shown 
in  Fig.  38.  The  velocity  rises  quite  sharply  from  its  initial  value  and 
reaches  a  maximum  at  about  2  kb  pressure.  As  pressure  is  further  in¬ 
creased,  the  velocity  slowly  decreases.  The  significance  of  this  feature 
is  not  clearly  understood,  and  its  validity  is  not  firmly  established. 

The  question  may  be  raised  that  errors  in  experiment  and  analysis 
are  great  enough  so  that  the  velocity  actually  reaches  s  constant  value 
or  perhaps  even  continues  to  increase  very  slowly,  as  is  the  expected 
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result.  However,  the  change  in  travel-time  plots  (Fig.  36)  indicates  that 
the  time  through  the  sample  and  pistons  reaches  a  nearly  constant  value 
after  passing  through  a  minimum.  The  time  through  the  pistons  is  ex¬ 
pected  either  to  remsin  constant  or  to  decrease  at  some  linear  rate. 

Also,  the  sample  length  decreases  linearly.  Therefore,  from  a  qualita¬ 
tive  point  of  view,  it  appears  that  the  maximum  in  the  compreasional 
veloci ty- stress  curve  (Fig.  38)  is  controlled  by  raw  travel-time  data 
and  not  by  errors  in  the  analysis. 

The  compressional  wave  travel-time  measurements  as  a  function  of 
pressure  are  thought  to  be  quite  accurate.  Changes  in  travel-time  can 
be  measured  to  ±0.005  /as  and  ram  pressure  and  lateral  pressure  can  be 
regulated  to  ±20  psi  The  uniaxial-strain  line  itself  is  not  well-known  but 
measurements  on  St.  Peter  sandstone  using  a  significantly  different 
uniaxial-strain  line  ahow  similar  results  for  change  in  travel-time. 

Sioux  quartzite,  Solenhofen  limestone,  Vermont  marble,  and  Massillon 
sandstone  do  not  show  such  a  minimum  in  the  travel-time  data.  For  other 


FIG.  38  ST.  PETER  SANDSTONE  -  COMPRESSIONAL  VELOCITY 
MEASUREMENTS 


M 


rocks,  time  continues  to  decrease  with  the  magnitude  of  decrease  becoming 
smaller  with  increasing  pressure. 

Shear  velocity  measurements  have  been  attempted  using  both  the  PSP 
technique  and  the  shear  transducers.  Reliable  measurements  as  a  function 
of  pressure  are  not  obtained  with  either  method.  The  most  serious  source 
of  trouble  is  that  an  increase  of  pressure  leads  to  a  severe  distortion  of 
the  waveform  being  observed.  The  PSP  measurements  are  especially  trouble¬ 
some  in  this  respect,  and  the  method  has  been  rejected  for  studying  porous 
sandstones . 

In  the  case  of  5  measurements  with  shear  transducers,  the  waveshape 
of  the  5  arrival  is  affected  by  changes  in  the  P  precursor  relative  to 
the  5  arrival.  The  P  precursor  is  very  small  at  its  beginning,  but  “noise" 
in  the  trace  increases  before  the  >S  arrival,  presumably  because  of  the 
arrival  of  waves  which  trsvel  part  of  the  distance  as  P  waves  and  part  as 
S  waves.  The  5  arrival  itself  also  appears  to  change  character  sometimes 
as  pressure  increases.  Getting  enough  energy  through  the  sample  ia 
another  source  of  difficulty  with  the  6'  wave  measurements — a  difficulty 
that  can  be  solved  with  a  stronger  pulse  generator.  The  first  problem 
will  still  remain  but  will  then  be  easier  to  overcome.  The  shear  velocity 
measurements  are  important  because  of  the  possibility  that  the  uniaxial- 
strain  ratio,  Tj j/Tj  j ,  may  change  with  pressure.  Such  a  change  may  be  the 
cause  of  the  compressional  velocity  decrease  with  pressure. 

2.  Massillon  Sandstone 

Four  runs  have  been  made  with  1-inch  specimens  of  Massillon  sandstone. 
Failure  occurs  at  4.2,  4.0,  3.4,  and  3.1  kb.  A  2-inch  specimen  failed  at 
4.0  kb.  The  uniaxial-strain  relationship,  t{2  ■  0.375  Tit,  which  had  been 
found  for  the  St.  Peter  sandstone,  has  been  used  for  these  measurements. 

It  is  thought  that  errors  involved  in  this  approximation  are  not  much 
greater  than  uncertainties  in  the  measurements,  but  for  future  work  the 
proper  relationship  should  be  developed.  Measurements  of  the  Hugoniot 
elastic  limit  give  values  of  6.9,  6.4,  and  4.2  kilobars,  with  the  lower 
value  being  favored  because  it  has  been  determined  from  a  record  of 
better  quality. 

The  dial  gauge  readings  are  less  linear  than  are  those  of  the 
St.  Peter  sandstone.  Some  evidence  of  yielding  has  been  noted  long 
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before  failure.  The  impurities  in  the  Maasillon  probably  make  the  stress- 
strain  relation  more  complex. 

The  congressional  wave  travel-time  behavior  appears  to  be  normal. 

The  minimum  that  characterized  the  St.  Peter  does  not  appear  in  the 
Massillon  except  possibly  just  before  failure.  A  sharp  initial  decrease 
is  gradually  attenuated  until  a  nearly  constant  travel-time  is  reached 
towards  the  end  of  the  run. 

3.  Coconino  Sandstone 

Two  runs  have  been  made  with  specimens  of  Coconino  sandstone  using 
the  St.  Peter  uniaxial-strain  relationship,  r2J  *  0.375  Tj j .  Both  failed 
at  about  3.9  kb.  Measurements  of  the  Hugoniot  elastic  limit  give  values 
of  'v2,  4.1,  5.1,  and  8  kb.  The  middle  two  values  are  more  reliable  be¬ 
cause  of  better  record  quality. 

Five  measurementa  taken  at  different  ratios  of  lateral  to  longitu¬ 
dinal  stress  provide  an  indication  of  the  error  involved  in  using  an 
incorrect  uniaxial-strain  relationship.  Using  the  ratio  TI2/Tn  "  0.23, 
the  specimens  fail  at  '*2. 6  kb;  and  using  the  ratio  0.47,  they  fail  at 
~4. 4  kb. 

The  dial  gauge  readings  are  more  linear  than  those  of  the  Maasillon, 
probably  reflecting  the  relative  lack  of  impurities  in  the  Coconino. 

(See  Section  IV- A- 2  for  more  detailed  description  of  composition.) 

4.  SOLENHOFEN  LIMESTONE 

Six  runs  have  been  made  with  two  specimens  of  Solenhofen  limestone. 
The  behavior  appears  quite  normal  with  the  compressional  velocity  in¬ 
creasing  from  5.56  mm/ixsec  at  zero-pressure  to  5.76  mm//xsec  at  5.2  kb. 

The  rate  of  increase,  as  expected,  lessens  with  pressure. 

Two  of  the  runs  have  been  made  with  shear  transducers.  The  behavior 
of  the  shear  velocity  is  similar  to  that  of  the  compressional  velocity, 
but  numerical  agreement  between  runs  is  poor. 

No  extensive  yielding  or  other  anomalous  behavior  of  the  rock  has 
been  noted.  This  is  in  contrast  with  the  Vermont  marble  which  displays 
plastic  yielding  in  the  range  3.0  kb  to  5.7  kb. 
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